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ANALYSIS FRAMES WITH NONLINEAR BEHAVIOR 


Alfredo Hua-Sing A.M. ASCE 


SYNOPSIS 


method for the analysis frames with nonlinear behavior subjected 
lateral deformations presented. The nonlinear characteristics both the 
members and the connections are considered. The relationships needed for 
the analysis mild steel frames including strain hardening are derived for 
use with any structural wide flange sections with error less than 
the resisting moments. Two examples are presented illustrate the applica- 
tions the method. 


INTRODUCTION 


Analyses indeterminate frame structures mild steel subjected de- 
formations beyond the elastic limit are usually based moment-curvature 
relationships that are flat after the fully-plastic resistance member has 
been developed. the deformations the members are not restricted 
buckling fracture, the behavior structure largely function the 
stress-strain characteristics the material, and strengths greater than the 
fully-plastic resistance can obtained result strain hardening. Often 
desirable able assess the benefits associated with such increased 
strength. example the case structure designed resist blast load- 
ing for which desired determine the energy-absorbing capacity when 
strain hardening taken into account. Likewise, research studies, 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Asst. Prof., Dept. Civ. Engrg., Univ. Urbana, 
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often desirable available means analyzing the resistance frames 
and beams that permits strain hardening considered. 

The subject elasto-plastic analysis and pure plastic analysis have been 
treated extensively several authors. However, very limited amount 
published material can found that treats the analysis statically indeter- 
minate structures when the nonlinear behavior its component parts con- 
sidered. This paper presents method which indeterminate frame struc- 
tures composed elements (members and connections) individual 
resistance-deformation characteristics any monotonically increasing form 
subjected lateral deformations can analyzed. For indeterminate struc- 
tures mild steel, the increased strength and energy-absorbing capacity re- 
sulting from strain hardening can, therefore, calculated. The resistance 
structures with riveted and bolted column-base and beam-to-column connec- 
tions can evaluated with the nonlinear characteristics the connections 
materials with nonlinear characteristics, such aluminum, without resorting 
linear approximations. The method has been found particularly useful 
the case research studies. 

The practicability the method depends upon means relating 
the resistance member particular deformation. this method the 
resisting moments the ends member are related the corresponding 
end-slopes. This relationship called the “moment end-slope relationship.” 

The resisting moments that correspond given lateral displacements 
the loaded joints structure are found means iterative procedure. 
Sucha procedure made convenient the use the curves 
for the individual members and moment-rotation curves for the nonrigid con- 
nections. After the compatible set resisting moments have been obtained, 
the corresponding set loads required produce the particular joint displace- 
ments are computed. solving set such problems, load-joint displace- 
ment relationships can obtained for range loads, conversely, for 
range displacements. 

particular, the application the method the analysis mild-steel 
frames loaded into the range strain hardening presented. Two examples 
are solved illustrate the application the method. The structure Example 
rigid joints, while that Example nonrigid joints. 

assumed that the ultimate strength mild steel direct tension 
compression test, under the particular time-dependent conditions involved, can 
fully utilized the structural members. The stress-strain relationship 
obtaining the moment end-slope relationship typical that obtained 
for ASTM A-7 steel, shown Fig. The moment end-slope relationship 
for mild steel are presented dimensionless form and can used for all 
commonly available structural wide flange sections within error 
the resisting moments. 

Notation.—The letter symbols adopted for use this paper are defined where 
they first appear, the illustrations the text, and are arranged alpha- 
betically, for convenience reference, the Appendix. 


THE METHOD ANALYSIS 


The assumption underlying the method analysis are follows: 
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The moment-curvature relationship corresponds pure flexure. Effects 
axial forces, shearing forces, buckling, and residual stresses such re- 
lationship are neglected. 

Loadings are concentrated the joints and are always increasing. 

Clockwise moments acting member are positive. 

joint, moments and including the ultimate moment 
the members can transmitted without destroying continuity the joint. For 
anonrigid joint, the moment-rotation characteristics the connection must 
defined and including the ultimate moment. such cases, the rotations 
both the connection and the members must considered. 


The analysis statically indeterminate structures involves, basically, the 
solutions the equations equilibrium andcontinuity. The equations equi- 
librium and physical continuity joint Fig. 2,in given structure may 
expressed, respectively, as: 


and denotes the number members meeting joint, the total moment 
section, and the total rotation. Ifa simple relationship between the mo- 


ments, and the end-slopes and connection rotations can found, the total ro- 


and equations continuity, conversely, equations equilibrium, may 
written for each joint structure that are analogous the slope-deflection 
equations for linearly elastic structures. Under such conditions purely ana- 
lytical solution possible. For nonlinear cases, the relationship defined 
Eq. usually turns out complicated, the case mild steel with 
strain hardening, and the equations equilibrium continuity will involve 
integralexpressions. For such cases, semi-graphical step-by-step approxi- 
are expressed graphical form. the end slope end member ji, 
and the rotation the connection member ji. 

The moments the ends every member structure are assumed 
random. These moments, which are general not compatible with given de- 
formation pattern the structure, are corrected iteratively, subject the 
conditions equilibrium and continuity the joints, until set values with 
convergence desired degree precision has been obtained. starts out 
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assuming displacements the loaded joints. From this assumption the 
general deflected shape the structure may determined, and the finding 
the appropriate resistance the structure that particular deformation pat- 
tern accomplished, assuming that reversal not occur any place 
the structure. This requires prior knowledge the correct direction 
the moment, determined from the deflected shape the structure. 

The iterative method correcting the end moments based the prin- 
ciple that the moment the near-end member, shown Fig. 
can found that will produce the the near-end with 
known assumed far-end moment, The near-end moment, that 
corresponds given end-slope, and far-end moment, can de- 
termined rapidly from the moment end-slope curves given Figs. and 
(Fig. portion Fig. enlarged scale). Figs. are for mem- 
bers with contraflexure; for members without contraflexure. The val- 
ues the connection rotation, may determined from moment-rotation 
curves Figs. Figs. and are presented herein for pur- 
poses illustration only, but are typical actual beam-to-column 
base moment-rotation relationships. 

The joint obtained once set near-end moments and the 
corresponding near-end rotations found such that the relationships defined 
Eqs. and are satisfied simultaneously. 

The moments and rotations joint would the correct values all the 
far-end moments, that were used determining and were the 
correct ones. The determination this tentative solution for joint suc- 
cessively carried out for all the joints given structure, with the most re- 
cently found moments the opposite ends the members the far-end 
moments for the joint concerned. The final solution for structure found 
when the tentative solutions ofevery jointare the correct solutions the joints. 
Convergence quite rapid and, ordinarily, not more than three cycles are re- 


quired obtain accuracy consistent with the accuracy the moment end- 
curves. 


MOMENT-CURVATURE RELATIONSHIP 


The moment-curvature relationship for mild steel has been determined 
dimensionless form such that the same relationship applies all commonly 
available structural wide flange sections A-7 steel with maximum possible 
error resulting from variation geometrical properties less than 3%in 
the resisting moments. 

The derivation the dimensionless moment-curvature relationships 
based the following assumptions: 


wide flange section can idealized being composed three rec- 
tangular components. 

Strains are distributed linearly through the depth section. 

The stress-strain relationship for the material the same tension 
and compression, and the resistance member flexure directly de- 
pendent upon this relationship. 


The average stress the flange, the stress the centroidal axis 
the flange. 
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FIG. SLOPE RELATIONSHIPS FOR SECTIONS ASTM STEEL 
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FIG, 6.—ASSUMED MOMENT-ROTATION 
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With these assumptions, the dimensionless moment may expressed 


terms the elastic limit moment section, 


which factor determined the yield stress level and stress distri- 

this form, study all commonly available structural wide flange sec- 
tions A-7 steel indicates maximum possible difference M/M, about 
12%. average moment-curvature relationship, shown can used 
for all commonly available wide flange sections A-7 steel with maximum 
possible error the resisting moments. 

The variation the values closely related with the variation 
the shape factors, which the ratio section. The shape factor, 
S.F., wide flange section given 


d-t 


The resisting moment cantherefore expressed terms the fully-plastic 


this form, study all commonly available structural wide flange sec- 
tions A-7 steel indicates maximum possible difference about 
The average moment- relationship, shown Fig. can used 


for all structural wide flange sections with maximum possible error less 
than the resisting moments. 


MOMENT END-SLOPE RELATIONSHIP 


The method analysis utilizes the moment end-slope relationships the 
members. These relationships are those simply-supported beam with 
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moments applied both ends. The derivation these relationships based 
the following assumptions: 


Anend slope positive the direction its rotation the same 
the end moment. 

Relative lateral displacement described the product angle- 
change times the original length member. 


With these assumptions, the dimensionless moment end-slope relationships 
the total angle change beam with applied one end, 
Fig. 9(c). 

The simply-supported beam with moments applied both ends shown 
Fig. 9(a). The moment diagram (in dash lines) and its corresponding curva- 
ture diagram (in bold lines) are shown Fig. 9(b). The curvature diagram 
deduced from the average moment-curvature relationship Fig. The end 
Slopes given set ofend moments are the corresponding reactions the 
supports with the curvature diagram acting the load the beam. This can 
expressed with reference Fig. 9(b) 


L-x 


and 


These end-slopes can obtained easily using numerical 
and considering the simply-supported beam Fig. 9(a) being cantilevered 
and free The total angle change between and and the perpen- 
dicular distance from the line tangent the deflected beam Ajj, are 
readily obtained numerical integration. The end slopes are, therefore, 


and 
was shown earlier that average moment-curvature relationship 
terms the fully plastic moments can used for all commonly available 
structural wide flange sections with maximum possible error less than 
the resisting moment. The moment-end slope relationships given 
Figs. and which were derived with the average moment-curvature curve 
Fig. can shownto have maximum possible error resulting from vari- 


ation geometrical properties less than relationships are 
used for all commonly available wide flange sections A-7 steel. 


Procedure for Computing Deflections, Moments, and Buckling Loads,” 
Newmark, Transactions, ASCE, Vol. 108, 1945. 
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9.—DESIGNATION MOMENTS AND DEFORMATIONS 
SIMPLE SUPPORTED BEAM 


| 
| 
1 — ~, - 2 
/ 
/ 
/ 
/ 
/ 
(a) 
~ 
~ 
| 
(b 
i 
a 
4 
} 
hes, 
ae 


The deflected configuration member, whether with without contra- 


ILLUSTRATIVE PROBLEMS 


The two sample problems illustrated the following sections are simple 
frames with rigid and nonrigid joints, respectively. However, the method ap- 
plies just well the analysis more complicated frames. Examples 
multi-story, multi-bay frames with rigid well multi-story, 
multi-bay frames with nonrigid joints3,5 may found elsewhere. 

Example Single Story, Single Bay Frame with Rigid Joints.—The struc- 
ture shown Fig. 10(a) subjected lateral deformations, ft. 
Determine the load corresponding these deformations. Properties the 
the rectangular moment inertia 2,942.4 367.8 Shape 
11,770k-in. 981k-ft. Mg, 11,770 1.16 13,653k-in. 1,138 k-ft. 

Since the members are 228 section, the fully-plastic moments 
are the same and will designated the calculations below. The 
elastic limit angle-changes for the members are: 


The chord displacement angle, shown Fig. 10(c) 


The deflected shape the structure and the correct directions the mo- 
ments are also shown Fig. 10(c). The initial assumedend moments are given 
the first line moments (nearest the member) Fig. 10(b) with clockwise 
moments the members considered positive. For the columns these mo- 
ments are the fixed end moments due the given displacements which are 
found from Fig. 3(a). For column AB, these are determined taking ji/ 
and thus giving value 1.61 

The first two lines the computations Table will explained detail. 


Referring Fig. 10(c), the equilibrium and continuity requirements joint 
are, respectively, 


Method for the Analysis Frames Subjected Inelastic Deformation into the 
Range Strain Hardening,” Ang, Thesis, Univ. Urbana, 
February, 1957. 

Method for the Analysis Frames Subjected Inelastic Deformation into the 
Range Strain Hardening,” Ang and Massard, Univ. Illinois Civ. Engrg. 
Studies, Struct. Research Series No. 130, November, 1956. 

“An Investigation Riveted and Bolted Column-Base and Beam-to-Column Con- 
nections Under Slow and Rapid Loading,” McDonald, Ang and Massard, 
Univ. Civ. Engrg. Studies, Struct. Research Series No. 156, February, 1958. 
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FIG, SINGLE STORY, SINGLE BAY FRAME WITH RIGID 
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Eqs. and may solved simultaneously trial-and- first 
trial assume the moments joint MBA MBC 1.20 Mfp These val- 
ues the moments automatically satisfy Eq. and will checked for Eq. 
follows. Since both members and have contraflexure, Fig. 3(a) 
3(b) should used for these members. For member BA, enter Fig. 3(a) 3(b) 
ment A); the end-slope corre these moments given Fig. 


the same figures with 1.20 Mg, and (the assumed 
moment member BC); the end-slope corresponding these moments 


rotation for member as, 8.00 The assumed mo- 
ments 1.20 Mf, satisfy Eq. but not satisfy Eq. since 
9.61 8.00 Therefore these values the moments are incorrect. 

From the magnitudes the total rotations given the inequality above, 
the same table show the second trial these same moments. With 


values the end moments Eq. satisfied and Eqs. are approximately 
satisfied since 9.11 9.13 Therefore, Mpa Mpc 1.23 are 
considered correct for this cycle andare second line moments 
joint Fig. 10(b). 

The complete moments and rotations given Table and 
the convergence moments shown Fig. 10(b), which the underlined 
values are the final end moments terms 

The load which equal the sum the shear forces the columns, 
may calculated 


1138 1.34) 1138 1.24 314.7 kips. 


Pec 

Example Single Story, Single Bay Frame with Non-Rigid Joints.—The 
same Example Awill solved considering that the joints are 
nonrigid (Fig. 11(a)) and have the moment-rotation relationships for beam- 
to-column and column-base connections given Figs. and respectively. 
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These moment-rotation relationships are assumed for purposes illustration 
only. For convenience, the moments are expressed, respectively, 
terms Mfp and the connected members. The columns are assumed 

The deflected shape the structure and the correct direction moments 
are shown Fig. 11(c). Initial end moments 1.00 are assumed for all 
members, except joint where the moment zero. These moments are 


TABLE 1.—COMPUTATION MOMENTS AND ROTATIONS FOR EXAMPLE 


CYCLE 


Ass. 1.23 
1.23 


Ass. 


CYCLE 


Cycle No.3 gives the same values those found Cycle No. Therefore the values 
Cycle are the correct moments. 


the values given closest the members Fig. 11(b). Clockwise moments 
acting the members are positive. 
The first two lines the computations Table will explained 


Referring Fig. 11(c), the equilibrium and continuity requirements joint 
are, respectively, 


(1) (2) (3) (4) (5) (6) 
Ass. 1.20 1.50 9.61 9.61 
1.20 6.00 6.00 8.00 
6.85 6.85 9.13 
Ass. 1.20 6.00 5.11 5.11 
1.20 2.90 2.90 3.86 
Ass. 1,23 6.75 4.36 4.36 
1.23 3.35 3.35 4.46 
Ass. 1.30 3.55 7.76 7.76 
1.30 4.55 4.55 6.06 
Ass. 1.34 4.10 7.01 7.01 
1.34 5.35 5.35 7.13 

3.15 3.15 4.20 
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FIG, 11.—SOLUTION SINGLE STORY, SINGLE BAY 
FRAME WITH NON-RIGID JOINTS. 
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TABLE 2.—COMPUTATION MOMENTS AND ROTATIONS FOR EXAMPLE 


CYCLE 


Ass. 1.30 5.10 2.90 
Ass. 1.40 7.20 4.00 
Ass. 
1.05 


CYCLE No. 


Ass. 1.06 
1.06 


Ass. 1.07 
1.07 


CYCLE No. 


Cycle No.3 gives the same values those found Cycle No. Therefore the values 
Cycle No, are the correct moments. 


Eqs. and may solved simultaneously trial-and-error. first 
trial assume the moments joint Mpa Mpc 1.05 For mem- 
ber BA, enter Fig. 3(a) 3(b) with Mpa 1.05 and Map 
1.00 Mfp (the assumed moment A); the end-slope corresponding these 


(1) (2) (3) (4) (5) (6) (7) 
Ass. 1.05 1.10 10.01 10.01 
1.05 1.10 5.80 6.90 9.20 
7.25 9.67 
9.01 9.01 
1.25 6.00 7.25 9.67 
1.40 6.30 7.70 10.27 
Ass. 1.05 1.00 5.70 6.70 


10.01 For member BC, enter the same figures with 
and Mcp 1.00 (the assumed moment for member BC). The 


end-slope corresponding these moments, given Fig. 3(a) 3(b), 
assumed moments Mpa 1.05 satisfy Eq. but not satisfy Eq. 


since 10.01 9.20 Therefore these assumed values the moments 
are incorrect. The correct values these moments are found 1.06 
for this cycle, shown lines and the computations Table 
The complete computation for moments and rotations given Table 
and the convergence moments shown Fig. 11(b) the underlined 
values are the final end moments given terms 
The load which equal the sum the shear forces the columns, 
may computed 


and the connection rotation for member 


267.0 kips 
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NOTATION 


The following notation has been adopted this paper: 


Superscripts: designates the member considered 
designates the connection considered 


designates quantity used the reference value 
problem 


Subscripts: designates elastic limit condition 
designates fully-plastic condition 

moment inertia about the centroidal axis the cross section 
overall depth structural section 

section modulus 

stress 

strain 

total moment section 

applied lateral load 

elastic limit resisting moment member 

fully-plastic resisting moment member 


fully-plastic resisting moment particular member used stan- 
dard unit moment problem 


total moment end member 

static yield stress the material 

total angle-change along the full length member 


total elastic limit angle-change member for moment applied 
one end. (See Fig. 9c.) 


total elastic limit angle-change particular member used 
standard unit rotation problem 


end slope end member 
chord displacement angle 


total rotation joing member 
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BEHAVIOR VISCOELASTIC PLATES BENDING 


George Masel 


SYNOPSIS 


Solutions the fundamental equation for linear viscoelastic plates bend- 
ing are obtained use the Laplace transformation. The quasi-static de- 
flection under various lateral loadings and the dynamic response under load 
are presented. Specific problems worked out include (a) the Kelvin and Max- 
well type plates under proportional loading, (b) rectangular Maxwell plate 


under moving load, and(c) free vibrations rectangular Maxwell and Kelvin 
type plates. 


INTRODUCTION 


number solutions viscoelastic problems involving tie concept 
differential viscoelastic operators have the 
Although these papers have been primarily concerned with stress analysis 
problems, the differential operator concept serves well problems for 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Prof. Applied Mechanics, Michigan State Univ., East Lansing, Mich. 


“Stress Analysis Viscoelastic Bodies,” E.H. Lee, Quarterly Applied Math- 
ematics, Vol. XIII, No. July, 1955, 183. 


“Stresses Elastically Reinforced, Viscoelastic Tubes with Internal Pressure,” 
Radok and Lee, Technical Report No. 15, Nord 16471, Brown Univ., 
April, 1956. 

Viscoelastic Anisotropic Media,” Biot, Proceedings, Se- 
cond Midwestern Conference Solid Mechanics, Purdue Univ., September, 1955, 94. 

“On Transient Thermal Stresses Linear Viscoelasticity,” Eli Sternberg, Pro- 
ceedings, Third National Congress Applied Mechanics, published 1958. 
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which deformations deflections are desired. this paper, the behavior 
flat viscoelastic plates bending studied. The basic differential equations 
are established using the operator concept and the solutions obtained ap- 
plication the Laplace transform. Two classes problems are considered. 
The quasi-static deflection laterally loaded plates obtained terms 
the associated elastic deflection similar plate. Plate moments for this 
case are also included and compared with those the elastic case. The 
second class considered that the dynamic response simply supported 
rectangular plates under load. Here, inertia forces due are 
taken into account and the free vibrations Maxwell and Kelvin type plate 
are given. Throughout the work physical constants the viscoelastic medium 
are assumed constant. 


FUNDAMENTAL EQUATIONS 


For isotropic medium linear viscoelastic behavior may defined 
means the linear differential equations 


here and designate the components the stress and strain tensors, re- 
spectively, and sjj and the components their deviators given 


which stands for the Kronecker delta. the prior equations and here- 
after throughout this paper, subscripts, unless otherwise specified, take 
the values and and the usual conventions for summation over repeated 
indices and for space differentiation are The operators and 
are linear differential time operators the form 


The coefficients etc. which describe the physical properties the 
material are assumed tobe temperature-independent and are, therefore, taken 
constants the following work. 

Experimental evidence indicates that most engineering materials behave 
elastically under hydrostatic pressure. This permits immediate simplifica- 
tion Eq. and the viscoelastic stress-strain relations may given 


4 


VISCOELASTIC PLATES 


which the usual bulk modulus elasticity theory. This form the 
basic equations will adopted here. 

Certain specializations Eq. interest from the point view ap- 
plications. When, for example, and the shear modulus, the equa- 
tions define elastic material. Using operators first order the equations 
may put into suitableforms which define thefrequently encountered Maxwell 
and Kelvin (Voigt) materials. Thus, 


define the Maxwell body characterized instantaneous elastic response fol- 
lowed viscous flow. Here the dot denotes differentiation with respect 
time and the coefficient viscosity. Likewise 


define aKelvin body characterized delayed elastic response. The phenome- 
nological behavior such viscoelastic media those defined previously 
frequently associated with the response mechanical models composed 
spring and dashpot elements arranged suitable fashion. Such models are 
discussed Alfrey6 and Gross.7 

The fact that the stress-strain equations the theory elasticity are in- 
corporated special case the basic stress-strain relations visco- 
elasticity theory has lead the development analogies relating the equa- 
tions those the other. Alfrey’s for the incompressible 
case has been extended cover compressible media Tsien.9 
Lee2 has established analogy for isotropic using the Laplace trans- 
form. More recently, has given correspondence principle for 
anisotropic media based upon the formal analogy the operational tensor 
viscoelastic theory and the elastic moduli elasticity theory. These various 
analogies may interpreted schemes whereby the elastic constants 
given elasticity equation are replaced differential time operators arrive 
analogous viscoelastic equation. For example, may seen from 
Eq. the shear modulus corresponds the operator P/Q. 

the particular problem hand the operational form the plate flexural 


rigidity required. Recalling that for plate thickness 


“Mechanical High Polymers,” Alfrey, Interscience Publishers, 

New York, Y., 1948. 
Structure the Theories Viscoelasticity,” Gross, Her- 
mann and Cie, Paris, 

“Non-homogeneous Stresses Viscoelastic Media,” Alfrey, Quarterly 
Applied Mathematics, Vol. May, 1944, pp. 113-119. 

Generalization Alfrey’s Theorem for Viscoelastic Media,” Tsien, 
arterly Applied Mathematics, Vol. 1950, 104. 


q 

= 

. 

Pe 

: 

4 

a 

3 

Bic 

i 

3 

‘ 


June, 


equivalent but more appropriate form for this work 


3 (3K + 4G 0 & 


The operational form based upon Eq. given 


For the purpose later reference convenient list here the specific 
forms this operator for the Maxwell and Kelvin materials. Thus, for 
Maxwell material 


(9) 


and for Kelvin material 


(10) 


the following, the plates considered are assumed tobe thin plates con- 
stant thickness Let the plane, referred set right-handed rec- 
tangular Cartesian coordinates, the middle plane the plate. The positive 
direction taken downward and the deflection the median plane this 
direction given distributed load f(x,y,t) per unit area the x,y 
plane acts the direction. the plate density. Accordingly, the plate 
equation for elastic plates when inertiaforcesdue deformation are taken in- 
account 


Replacing this the operator gives the fundamental equation 
for viscoelastic plate response, namely, 


QUASI-STATIC SOLUTIONS 


Proportional inertia forces due deformation are neglected, 
Eq. reduces the basic equation governing the quasi-static deflection 
viscoelastic plates. This equation may written 


3 
4 
(8) 
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VISCOELASTIC PLATES 


When the load function, such that may given the product space 
function multiplied time function, the variables may separated and 
solution terms the associated elastic solution obtained. Thus, letting 


yields upon substitution into Eq. the equation 


this equation, the right hand side function and and the left hand 
side function therefore, each side must equal the same con- 
stant. Taking this constant the flexural rigidity the plate for pur- 
poses convenience leads the equations 


Eq. simply the elastic plate equation for which solutions are already 
known. Such solutions may expressed the form 


(x,y) 
Solutions Eq. may obtained conveniently with the aid the Laplace 
transform method the operational calculus, assumed that the plate 
rest its undeflected position and free from all stress the instant 
application the load, the operator will transform into rational function 
the transform variable. Therefore, the transform Eq. becomes 


which barred quantities indicate the transforms the same quantity with- 
out the bar and g(s) and k(s) are polynomials. may solved for and 


the function expanded partial fractions give 


i 


- 
= 
4 
(21) 
i=1 


which the roots are all Their number, will depend upon the 
degree the polynomial, Eq. may inverted once the convo- 
lution integral give the general solution Eq. 


-r.T 


which the Dirac delta function and the base the natural log- 
arithm system. From this and Eq. follows that 


particular example the previous solution, consider plate made 
Kelvin material and having load distribution for which the elastic solution 
known. the load assumed tobe applied suddenly time and main- 
tained constant thereafter, given the Heaviside step function defined 


The constants and roots may obtained from the transform 
given Eq. 10, and Eq. then integrated directly give 


3K+G 


the load the plate inthe previous example suddenly removed some 
time the expression for the deflection for given 


(25) 


Fig. shows plot w/we versus obtained from Eqs. and fora 
plate having 0.25 and for which 2T. 

similar manner expressions for the viscoelastic plate moments 
may obtained from the corresponding equations for the 
elastic plate moments given 


“Theory Stress-Strain Relations Anisotropic Viscoelasticity and Relaxation 


Phenomena,” Biot, Journal Applied Physics, Vol. 25, No. 11, November, 
1954, 1385. 
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For aKelvin plate such was used the prior example the shear modulus 


which may inverted once give 


Examination Eq. reveals that othe viscoelastic plate moments ap- 
proach the elastic ones and when the moments are those elastic 
plate having which corresponds the case rigid body. 

The analysis Maxwell-type plate under the same loading conditions 
would follow the same fashion that for the Kelvin plate previously 
given. Omitting the details, the equations are Maxwell plate corresponding 
Eqs. 24, 25, and are given 


and 
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Eqs. and for the Maxwell plate are plotted Fig. using the same con- 


ditions 0.25 and 27) employed plotting the Kelvin plate curves 
shown there. 


LOAD 


~ 
— 


FIG, 1.—MAXWELL AND KELVIN PLATE DEFLECTION 


Moving Line Load.—If the load, f(x, cannot expressed simply 
the product space function multiplied time function, the question 
separating the variables the fundamental equation not straight forward 
operation and modified approach required. Consider, for example, the 
problem line load constant magnitude moving with constant velocity 
across simply supported rectangular plate having the dimensions shown 
Fig. view the boundary conditions the edge the plate, con- 
venient take the deflection the form the series 


m=1,3,5 n=1,3,5 


The procedure now expand the load function, similar series and 
substituting these series into the fundamental plate equation arrive equa- 


tion from which the coefficients, the deflection expression ob- 
tained. this end let 


m=1 n=1 


e 
\ 
~ 
= 
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Considering the line load uniform intensity over the small interval 
the plate and assuming the load the left edge the plate and moving 


FIG, 2.—RECTANGULAR PLATE FIG, 3.—MOVING LOAD REPRESENTATION 


Eq. may determined the usual manner based upon the orthogonality 
the sine functions. For the case hand the relationship given 


vt-e 


Carrying out the indicated integrations and inserting the limits yields the re- 
sult 


m=1,3,5 n=1,3,5 


Substitution this series together with the one given Eq. into Eq. 
ind matching coefficients produces anordinary differential equation from which 


unknown Eq. may determined. The equation tobe 
iolved 
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from which operational form given 


Eq. may solved once using the Laplace transform method the 
previous work. If, for example, the plate assumed made Maxwell 
type material, the operator given Eq.9 may introduced into Eq. 
and the solution obtained direct application the operational procedure. 
Omitting the details, the result for such plate given 


undeflected plate before the load has moved upon it. 

essential note that the solution given Eq. valid only while 
the load actually the plate, that is, when <a/v. the instant the 
load leaves the plate the deflection given 


i 


VISCOELASTIC PLATES 
and the plate velocity given 


(41) 


Using these expressions initial conditions the deflection the plate subse- 
quent the passing the load may bedetermined the same method was 
used for the loaded plate. The operational equation solved here 


which for the Maxwell case has the solution 


which 


mn] 


t=a/v 

Finally, the deflection for a/v given 


m=1,3,5 n=1,3,5 
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would expected from the model representation Maxwell material 
the deflection does not reduce zero increases indefinitely, but some 
residual deflection remains the plate. This given 


DYNAMIC SOLUTIONS 


Compressible Plate Material.—If inertia effects due deformation are 
taken into account, the second term the left hand side Eq. must 


retained. the present work only free vibrations are studied. Accordingly, 
solutions 


are sought. 

Solutions Eq. may obtained using the method the previous sec- 
tion. Consider example the free vibrations simply supported rec- 
tangular plate the dimensions given Fig.2. Taking again the deflection 
the plate given the double sine series Eq. and substituting this 
into Eq. results the following equation for determining mn. 


the plate assumed Maxwell plate this equation becomes 


(48) 


which 
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VISCOELASTIC PLATES 
Taking the Laplace transform Eq. results 


which polynomial the transform variable, the exact form 
which depends upon the initial conditions the plate motion. Eq. may 
inverted the method partial fractions once the roots the cubic factor 
the denominator have been found. Standard procedures for determining 
these roots are available.11 


Calling the roots rj, and the solution Eq. becomes 


1 


which the are constants involving the physical parameters the 
plate material and initial conditions the motion. Hence, the response 
given 


i=1 


m=1,3,5 n=1,3,5 


From this, numerical values are known and the initial con- 
ditions are given, the plate response may expressed terms the normal 
modes. should pointed out that because the presence the 
cubic Eq. the roots will bedifferent for each mode vibration. There- 
fore, considerable calculation would encountered carrying out the expan- 
sion Eq. higher modes. 

Incompressible Plate the plate material may assumed 
incompressible, frequently done viscoelastic problems, considerable 
simplification Eq. 48is achieved and general solutions are readily obtained. 
This essentially the same considering Under this assumption the 
form the flexural rigidity becomes 


and Eq. now 


“Theory Equations,” Dickson, John Wiley, and Sons, Inc., New York, 
1939, 42. 
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Solutions Eq. may obtained once. Consider, for example, the case 
simply supported Maxwell plate. There loss generality assum- 
ing initial conditions such that plate has velocity but displace- 
ment acceleration. For this condition the transform Eq. given 


verting gives 


For this, that is, the “damping effect” the viscous nature 
the plate small enough the response will oscillatory. Hence, inverting 
the equation for this case gives the result 


mn 4 
which 
T 


Thus, the response consists constant term plus decaying sine term. 

equal interest the previous example from practical point view 
the free response Kelvin plate, which corresponds essentially the con- 
cept internal viscousdamping used often studying damped vibrations 
beams and bars. Assuming that such plate started vibrating giving 
initial displacement and releasing without velocity the expression for 


the solution becomes 


the 
‘ 
| 
| 
j 
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which 


should noted this solution that the higher modes will rapidly 
damped out and, therefore, the fundamental mode that most impor- 
tance. 
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AIRCRAFT STRUCTURAL ANALYSIS ANALOG COMPUTER? 


SYNOPSIS 


has been many elastic structures can represented 
electric analogs with passive elements. recent years this methodof analy- 
sis has been successfully applied complex, highly redundant aircraft struc- 
tures. Presented this paper brief this technique andsome 
examples its practical application structural analysis. 


INTRODUCTION 


The structural design refinements required modern aircraft make 
necessary tohave and stress data under staticand dynamic 
loads. Present and future design trends require the analysis highly redun- 
dant structures. The principal analytical means use handle these prob- 
lems are either the formulating and solving largenumber simultaneous 
equations digital computer the use electric analog computers the 
direct analogy type. Digital computers have experienced rapid development 
recent years and solution time and solution accuracy are longer impedi- 
ments large scale structural application. The capabilities the electric 
analog method may presently not well established. The purpose this 
paper describe the use anelectric analog computer structural prob- 
lems Convair-Fort Worthand present the some applications 


Note.—Discussion open until November 1960. extend the one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, 86, No. June, 1960. 

Presented the October 1959 ASCE Convention Washington, 

Proj. Structures Engr., Convair, Div. General Dynamics, Fort Worth, Texas. 

Senior Structures Engr., Convair, Div. General Dynamics, Fort Worth, Texas. 
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order show that very satisfactory solutions extremely complex problems 
can obtained with this computer. 


GENERAL DESCRIPTION 


direct analogy electric analog computer may defined collection 
electrical elements (resistors, capacitors, transformers, etc.) that are used 
form circuits analogous the structure analyzed. Solutions are then 
obtained testing, taking measurements on, the electrical model much 
the same way that scale model would tested. Its principal features are: 


(1) furnishes the user with readily changeable model and is, therefore, 
most useful optimization and parameter variation studies. 

(2) like erector set that the elements can used over and over 
again different arrangements represent different structures. 


The use analogous electrical meansof solving pin-connected 
and rigidly-connected structures was proposed3 Bush 1934. Several 
other analogous circuits for isolated problems also appeared 
the literature prior 1945. that time the rather limited interest this 
method analysis was due the belief that satisfactory answers could ex- 
pected only for relatively simple problems because the inherent imperfec- 
tions the basic passive computer elements. Extensive efforts recent 
years the California Institute Technology, Pasadena, Calif., and Com- 
puter Engineering Associates Pasadena, Calif., have resulted the design 
computer that permits accurate solution large-scale problems. This 
computer design described This reference also lists large 
number papers the many analogies which have been developed for appli- 
cation this type computer. 

The computer usedin this work, designed for both static and 
ses, illustrated Fig. The computer contains the following: control 
desks; 154 inductors; 320 resistors; steady state current generators; 205 
transformers; general purpose amplifiers; and 164 date the 
computer has been used, alternately, perform: 


(1) Static stress analyses aircraft wings, fuselages, and tail surfaces. 

(2) Static aeroelastic analyses wings (loads are functions deflections). 

(3) Dynamic analyses including normal mode analyses, flutter gust re- 
sponse analyses entire aircraft and individual surfaces. 


This report will limited discussion the use the computer the 
static stress analyses field. 


BASIC CONCEPTS 


The Nodal Analogy.—The force-current, nodal analogy, normally used 
when solving structural problems the electric analog computer. the most 
general application dynamic problems, force analogous current and 


Analysis Electric Circuit Analogies,” Bush, Journal the 
Franklin Inst., March, 1934. 

Analogy for Shear Lag Problems,” Newton, Experimental 
Stress Analysis, Vol. II, No. October, 1944. 

“The Direct Analogy Analog,” McCann, ISA Journal, April, 1956. 
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voltage analogous velocity. For the case purely static stress prob- 
lems, analogous circuits using only resistors and transformers, which volt- 


age analogous displacement, are used. The basic relationships simple 
spring tension are summarized the Fig. 


Exact relations are dependent upon the choice scale factors. The choice 
the following basic factors may used. 


which and are constants, denotes electrical power and voltage. 
the general case multi-degree freedom system, the all 


DYNAMICS 
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FIG, 1.—LAYOUT CONVAIR DAEAC 


parts whereas may different for every coordinate. the current and 
the resistance. The choice form for the scale factor largely con- 


venience structural problems and the following derived expressions 
for the example Fig. 


The constants and are chosen one the last steps preparing 
problem for solution. They are chosen give the most practical size cir- 
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FIG, 2.—SPRING TENSION 


FIG, 3.—BAR TENSION 


(a) 


FIG, ELEMENT AND ITS ELECTRIC ANALOG 


cuit elements and also minimize parasitic effects. more complete dis- 
cussion scale factors may found eisewhere.® 

One-Dimensional Members continuous one-dimensional 
bers (straight curved bar tension, torsion, bending) the only approxi- 
mation necessary that distributed loads must replaced “equivalent” 
concentrated loads. For example, the spring previously illustrated represents 
exactly rod tension the value the spring constant, shown 
Fig. 

additional example one-dimensional member straight beam 
bending. Many beam analogies varying degrees accuracy are general 
use. The analogy follow exact under the limitations previously applied 
loads and equivalent the one presented Russell and Mac- 
although the derivation different. 

Shown Fig. 4(a) unloaded segment beam length, with in- 
ternal loads and coordinates defined the ends. The motion the right end 
relative the left end may expressed 


which Aj, Ag, are the flexibility influence coefficients for the right end 
the beam with the left end assumed fixed. indicated Fig. 4(a) the 
vertical force and the moment. 

The equilibrium equations are 


Now consider the circuit Fig. 4(b). Circuit equations are 


“The Solution Aeroelastic Problems Means Electrical Analogies,” 
MacNeal, McCann and Wilts, Journal Aero Sciences, Vol. 18, Decem- 
ber, 1951. 

Improved Electrical Analogy for the Analysis Beams Bending,” 
Russell and MacNeal, Journal Applied Mechanics, September, 1953. 
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Standard scale factors used the analysis are 

and 


The “P” factor used Eqs. and equivalent tothe use second con- 


stant Substituting Eqs. through into either set preceding 
equations gives 


a 
(18) 
and 


represent entire beam, common coordinates adjacent segments are 
joined. Supports and loads are introduced into circuits the appropriate co- 
ordinates. The “exactness” the analogy thus formed dependent upon the 
accuracy with which the influence coefficients can evaluated. noted 
that, the case beam with strain energy due only bending, the coef- 
ficients may evaluated from 


(22) 


which the second moment one area about right end. the 
the unit moment diagram, the deflection due the unit 
moment, and the deflection due the unit force. Additional items such 
shear deformationor the effect tapered flanges may handled includ- 
ing these effects the calculation the influence coefficients. 
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Two-Dimensional Elements —Aircraft structures are usually mixture 
one-dimensional elements (spar caps, stiffeners, fuselage frames rings) and 
two-dimensional elements such the flat curved web and cover panels. 
The treatment the two dimensional elements requires additional approxi- 
mations over the substitution concentrated loads for distributed loads. These 
approximations involve the lumping distributed stiffness elasticity into 
effective one-dimensional elements, and the substitution finite number 
interconnections for the continuous interconnections boundaries. The result 
that the actual structure replaced idealized one consisting ele- 
ments which individually can simulated. The final step the interconnection 

the analogies for the individual elements. 

The approach previously described, knownas the “lumped parameter” meth- 
od, generally more satisfactory structural problems than alternate 
“finite difference” method. this second approach the differential equations 
are first replaced finite difference equations, and circuits analogous the 
difference equations are derived. 

The treatment two-dimensional elements may illustrated flat 
plate plane stress. The plate shown Fig. replaced rectangular 
grid equivalent bars which carry axial load, and panels between the bars 
which carry shear loads only. Further, the panels are assumed connect- 
the bar eachedge the single point connection shown. Free body dia- 
grams the bars intersecting point are shown Fig. 6(a) and panel 
FLQK are shown Fig. 6(b). 

Consider first the bars intersecting For the assumed state 
uniform stress the neighborhood point the stress-strain relations are: 


Converting strain total elongation and writing stress terms forces, the 
equations become: 


A Uy Uy, €y Ay = Et (- Vv F Fy Ax ) ee ee . (27) 


The analogous circuit shown Fig. 7(a). The requirements for the circuit 
elements may established writing the circuit equations 
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(a) (b) 


6.—FREE BODY DIAGRAM BARS INTERSECTING 
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FIG, PLATE 
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Substitution the scale factor relations listed into Eqs. and 29, 


and letting gives the following, comparison with the mechani- 
cal Eqs. and 27: 


Analogous circuits for the shear panels will developed next using energy 
methods. During the deformation the panel, the external work done the 


edge forces must equal the strain energy stored the panel. The external 
work 


Workexyt = Fxy + Fyx (Oy, - Uy, ] 


The strain energy also can expressed terms 


Equating 


Fig. 7(b) contains the analogous circuit. The turns ratio satisfies the equi- 
librium equation: 


Substituting scale factors 

This compared with the mechanical equation 


2 
AE = I - 
Eta 


from which 


The electrical equation corresponding Eq. 


Substitution scale factors and the value from Eq. gives 


The combined circuit diagram for portionof the plate shownin Fig. given 
Fig. sections five bars and thetwo center shear panels. 

The process replacing the actual structure idealized one based 
mostly judgment. There nosingle method which applicable allstruc- 
tures. Each must treatedon individual basis and many factors including, 
large problems, the size the computer must considered. important 


FIG, 8.—FLAT PLATE ANALOGY 


feature the analog computer that reasonable representations general 
two and three-dimensional structures can built methods similar 
those used the preceding example. has not been the intent this paper 
compile extensive list usable analogies. The main purpose has been 
illustrate fundamentals which are related the large scale analyses 
discussed next. The interested reader referred other for 
comprehensive coverage analogies. 


“Electrical Analogs Statically Loaded Structures,” F.C. Ryder, Transactions, 
ASCE, Vol. 119, 1954. 
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ANALOG COMPUTER 
EXAMPLE BUILDING FRAME PROBLEM 


The frame structure shown Fig. will illustrate the use the analog 
computer composite structures composed one-dimensional elements. 
This example will also serve demonstrate the numerical accuracy which 
possible comparison with the solution obtained large 
scale digital computer. 

The frame Fig. will represented the proper connection six 
beam segments the type previously discussed. The connection the beam 
segments can done inspection the structure follows: 


(1) The slope beam point identical the slope beam 
point 

(2) The horizontal deflection beam point identical the hori- 
zontal deflection beam point and forth throughout the frame. 


FIG, 9.—EXAMPLE BUILDING FRAME PROBLEM 


These identities slopes and deflections may accomplished the analogy 
connecting the corresponding points the beam segment analo- 
gies. can easily shown that static equilibrium the joints the frame 
also maintained the analogy the connections discussed above. 

The circuit Fig. shows the electric analogy the frame Fig. 
The beam segment analogies have been kept separated that they may in- 
dividually observed. should noted that the number transformers re- 
quired the circuit Fig. could reduced six combining some 
the transformers shown. This reduction equipment would way effect 
the solution. 

The transformer ratios the circuit Fig. may combining 
Eqs. 18, 20, and 22. For constant 


Ty 


“Digital Computations for Matrix Analysis,” Archer, Proceed- 
ings, ASCE, Vol. 84, No. October, 1958. 
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The values the resistors Fig. may found combining 17, 
19, 21, and 23. 


Rj El (44) 


The numerical values the transformer ratios and resistors the circuit 
Fig. are evaluated Tables and The scale factors used are: 


C=1.5 104 
100 


The results the theframe are shownin Table along 
with the corresponding results obtained from Archer.9 


TWO-DIMENSIONAL STRUCTURAL ANALYSIS 


two-dimensional analysis plate-like structures, made vertical 
webs with cover skins the upper and lower surface, the assumption made 
that the structure symmetrical about central, chord plane, and that the 


BEAM BEAM G-D BEAM D-E BEAM E-F 


FIG, 10.—ANALOG CIRCUIT-BUILDING FRAME PROBLEM 


surface stresses and deformations are antisymmetrical about the same central 
plane. Forces and deformations the cover skins can then represented 
moments and rotations, respectively, about the central plane. The result 
that the skin representation reduced systemof beams and torque boxes 
acting conjunction with second beam network, representing the spars and 
ribs. 
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The analysis triangular airplane wing the analog computer that was 
used this work goodexample large scale analysis this type. The 
planform arrangement the vertical webs the structure shownin Fig. 11. 
The structure covered stiffened skin panels the top and bottom. 

Due the excessive size the overall problem, the structure Fig. 
was reduced complexity the one shown Fig. for initial analysis. 
The cover skins the structure were rectangular grid equiva- 
lent barsand shear discussed the section two-dimensional ele- 


ments. The shear panels are connected the spars and ribs the points 
shown Fig. 12. 


TABLE CALCULATIONS 


Length, Resistance Resistance 


TABLE 2.—TRANSFORMER CALCULATIONS 


Length, Number Number 
inches 
turns turns 


may noted Fig. that many these connections are not the 
periphery shear panel. For these cases modification, which considers 
the structural element shown Fig. 13, made the conventional shear 
panel analogy described under two-dimensional elements. Fig. 13(a) indicates 
skewed spans with orthogonal skin lumping and Fig. 13(b) indicates shear 
panel with interior forces. 

order obtain more detailed stress information the portion the 
structure Fig. aft line second analysis was made that portion 
the structure only. The used this secondanalysis are 
identical nature those the first. The structure was, however, repre- 
sented much finer detail due the reduction overall size. Fig. shows 


A 
120 400 1000 120 
120 400 1000 120 
240 400 2000 960 
240 320 2500 1200 
bay 
Beam 


the structural representation used this analysis. Along the line BB, where 


the wing “cut” for this analysis, boundary conditions obtained from the first 
analysis were applied. 


Extensive comparisons results from the second analysis with experi- 
mental results for several load conditions show good agreement. Analysis and 


TABLE FROM BUILDING FRAME ANALYSIS 


Analog Solution Digital Solution 
Point Vertical Horizontal Slope, Vertical Horizontal Slope, 
Deflection, Deflection, Deflection, Deflection, 
5.120 
5.120 
5.120 .0303 
1.91 5.120 


TABLE LOADS FROM BUILDING FRAME ANALYSIS 


Analog Solution Digital Solution 


Moment Moment Moment Moment 


experimental test results for typical condition are shown Figs. and 16. 


Items presented for comparison are skin span-wise normal stress and spar 
cap stress. 


THREE DIMENSIONAL STRUCTURAL ANALYSIS 


The analysis the inboard portion swept wing, also performed this 
analog computer, demonstrates the utility this type computer large 
scale three-dimensional structure. The idealized representation the struc- 
ture analyzed shown Fig. 17. 

The structure includes the inboard portionof the airplane wing and simpli- 
fied representation portion the fuselage. The fuselage structure was 


not the subject this analysis, but was included this analysis give proper 
support the wing. 
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FIG, 11.—TRIANGULAR WING 
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Beam with Web 
Beam without Web 


Shear Panel Bounded 


Skin Flanges 


Point Connection 


between Skin and 
Vertical Web 


FIG, 12.—REPRESENTATION FOR FIRST WING ANALYSIS 
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The wing structure consists three mainspars and number ribs which 
intersect non-orthogonally and have considerable amount taper. Further, 
the cover skins are curved and the structure lacks any degree symmetry 
about chord plane. Accuracy analysis results required individual repre- 
sentation all principal structural elements. Idealization consisted first 
replacing the cover skins number flat planes which approximated the 
curved surfaces. The treatment the cover skins which seemed most suit- 
able was, this instance, lump the skin tensile properties into bars fitting 


LUMPED 
SKIN BAR 


SHEAR PANEL 
ATTACHMENT 


POINTS RIB CAP AND 


LUMPED 


SPAR CAP 


(b) 
FIG, STRUCTURAL ELEMENT 


into the non-orthogonal spar-rib pattern. The bars which coincide with either 
spar rib cap were then lumped with the cap member. The effective bar 
pattern for the wing, along with the effective longerons the fuselage, are 
given the solid line pattern Between these bars lie the skewed 
cover skin shear panels and the shear panels representing the webs the spars 
and ribs. Typical skin and web panels are shown Figs. 18(a), skewed skin 
element, and 18(b), spar rib element. 

Because the span-wise and chord-wise axially loaded bars are not ortho- 
gonal each other, was considered necessary account for the coupling 
which exists between axial stresses and these bars and the shearing 
stress and strain the cover panels. This coupling can best described 
considering the elemental area skin shown Fig. 19. 
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FOR SECOND WING ANALYSIS 
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FIG, 18.—TYPICAL SKIN AND WEB PANELS 
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FIG, 17.—STRUCTURAL REPRESENTATION FOR SWEPT WING ANALYSIS 
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FIG, 21.—COMPARISON SPANWISE LOWER SURFACE STRESS 


The stress-strainrelationships Fig.19 are given Eq. 47, which shows 
that the shear and are not independent rec- 
tangular element. 


and are axial strains and denotes the shear strain. 

The analogy used represent the cover skins this analysis derived 
from theforegoing equation. The effect the coupling amounts 
connection between the panel and the surrounding bars, and illustrated 
the skewed element Fig. 18(a). 

The purpose this analysis was determine the structural deflections and 
internal loads due number different loading conditions. Each loading 
condition included calculated internal loads the outboard wing section this 
analysis. secondpurpose this analysis was the internal loads 
the structure due given loading condition after the presumed failure 
any chosen single structural member. This portion the analysis used 
demonstrate the safety the structure the event unobserved fatigue 
failure. The failure one structural member and quickly simu- 
lated this analysis effectively removing the portion the circuit repre- 
senting that member. 

Limited comparisons the results the analysis with available test data 
show very good agreement. Figs. and show comparison some nor- 


mal stresses the upper and lower cover skin areas ABCD and A'B'C'D' 
Fig. 17. 


The direct analogy electric analog computer has been found very 
useful toolin the analysis highly redundant structures. emphasized that 
the authors’ work the analog computer and the high speed digital computers 
are both extensively used for structural analysis. Each has class prob- 
lems for which best suited. The analog computer has been used princi- 
pally for the following classes problems: 


(1) Ondesign problems when many changes refinements may required 
arriving optimum arrangement. 

(2) Onsome problems for which satisfactory digital solution was not readi- 
available. these cases the use the analog computer has resultedin sav- 
ings both time and costs. 
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VIBRATIONS AND STABILITY PLATES UNDER INITIAL STRESS 
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SYNOPSIS 


theory the motion elastic plates established which, ad- 
dition initial membrane forces, takes into account initial moments and 
transverse shear forces. The equations motion are obtained linearizing 
the appropriate nonlinear plate equations and, the absence initial stresses, 
they are reduced, depending the order approximation, either the clas- 
sical(elementary) the Mindlin plate theory (classical theory extended in- 
clude the effects transverse sheardeformation and rotatory inertia). These 
equations are applied the particular problems and the effect 
uniform initial compression, bending moment and transverse shear force 
the frequencies the flexural, longitudinal, thickness-shear and the first two 
face-shear modes discussed. 

the course this investigation was found that plate under uniform 


transverse shear force may buckle fashion similar plate under com- 
pression. 


INTRODUCTION 


well-known that state initial stress may significantly change the 
static and dynamic response elastic body. One the most familiar 
amples concerns flexural vibrations hinged elastic beam. The funda- 
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mental frequency free vibrations diminished initial compression 
and may even vanish, the initial compression reaches the value the buck- 
ling stress. 

General, three-dimensional theories elasticity body under initial 
stress were formulated Biot (1)3 and Goodier (2). variety 
other studies have dealt with particular problems. For example, the effect 
initial stress the torsional stiffness was considered Biot (3) and 
Goodier (2), and the influence initial stress elastic waves Biot (4). 

plate theory, the interest the influence initial stress long stand- 
ing. The equilibrium plates, under initial stresses uniformly distributed 
across the plate thickness, has been studied Saint-Venant (5, 6), who 
derived the governing equations considering isolated plate element. 
correlation between Saint-Venant’s plate equations and the three-dimensional 
theory has been established (7). This existing theory plates limited 
initial “membrane” stresses. The possibility initial moments initial 
transverse shear forces has not, yet, been considered. 

appears that more complete theory equilibrium and motion plates, 
subjected more general state initial stress may useful, for should 
describe phenomena not contained Saint-Venant’s theory and should provide 
additional insight tothe reasons why only initial “membrane” stresses were 
considered previous investigations. 

addition generalizing the Saint-Venant plate theory, the intent this 
study establish more inclusive theory plates underinitial stress, an- 
alogous the Mindlin plate theory (8), retaining the equations 
motion the effects transverse shear deformation and rotatory inertia. 

Notation.—The letter symbols use inthis paper aredefined where 
they first appear, the illustrations the text, and are arranged alpha- 
betically, for the convenience reference, the Appendix. 


COMPONENTS STRAIN 


The equations motion for the isotropic elastic plate, that include the in- 
fluence initial stresses, will obtained linearizing non-linear plate 
equations, that will turn derived from basic relations non-linear elas- 
ticity. The form these non-linear equations dependent upon the specifi- 
cation stress strain and upon the order approximation. With regard 
strain, for example, Biot (1) has established second order theory wherein the 
strain components are linearly related the actual change length, rather 
than the change the square the length—as are the classical Lagrangian 
Eulerian components strain (see, for example, Novozhilov (9) and 
Sokolnikoff(10)). The results obtained using Biot’s theory are not sim- 
pler, however, than those obtained using the classical non-linear theory 
elasticity with Lagrangian strain components (15). For this reason, the latter 
will utilized this work. 

define the components strain, consider material element whose 
Cartesian coordinates before deformation are and whose coordinates 
after deformation, the same Cartesian system, are The components 


Numerals parentheses, thus (1), refer corresponding items the Bibliography. 
(see Appendix 
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displacement uy, are also dependent time, and are defined 
=z + uz (x, y, Zz, t) 


and the Lagrangian components strain 


DEFORMATIONS 


nate system, with the plane being the middle plane the plate. 

First, theory plates will established that includes the influence 
initial stresses, and that could reduced Mindlin’s theory (8) (classical 
theory, extended include the effects transverse shear deformation and 
rotatory inertia) the initial stresses vanish. 

has been done many previous investigations, the displacement com- 
ponents ux, Uy, are expanded infinite series powers the thickness 
coordinate and approximate displacement components Uy, Uz, are intro- 
duced truncating the infinite series. Only those terms are retained that 
contribute the modes vibration that will investigated. 

For this reason will suffice consider the following expansions; 
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STRAIN ENERGY 


consider abody original volume some general deformed position, 
specified displacement components ux, and shall calculate the incre- 
ment Wint the strain energy due arbitrary increments 
the displacement components. Since the strain energy function 
the strains, (more specifically the strain invariants) have 


energy density with respect the corresponding components strain 


Since the significance not immediately apparent, necessary 


consider the geometry deformation. general, three initially perpendicu- 
lar line elements dx, dy, dz, are displaced and strained due the deformation 


FIG, 1.—DEFORMED ELEMENT AND COMPONENTS STRESS 


into non-orthogonal parallelepiped. The force vector any one face this 
parallelepiped may resolved into (non-orthogonal) components parallel 
its edges, which, divided the face area before deformation, define the stress 
components obvious that the components not form symmet- 


ric stress tensor, since the moment condition equilibrium (Fig. required 
that 


the relative elongation line element that, before deformation, was 
parallel the x-axis. 
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symmetric stress tensor may formed the relations 


can shown (9) that which are Trefftz’s components stress (2), are 


Using the approximate expressions (Eq. 3), approximate increments the 
components strain (Eq. are established and substituted into Eq. Since 
the dependence the z-coordinate specified, integration through the thick- 
ness, readily carried out. This process leads the plate stresses 


h/2 +h/2 
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The surface integral extends over the face area the plate before de- 
formation; the line integral extends over the total undeformed cylindrical 
boundary the plate; and are coordinates measured along the outward 
normal and along the tangent the undeformed boundary, respectively. Notice 


that Eq. the contribution the plate stresses Nzz, Kxx, Tx, Ty, 
the increment the strain energy, disregarded. 


WORK EXTERNAL FORCES 


The external forces consist boundary tractions, taken for convenience 
per unit undeformed boundary area, and body forces taken per unit unde- 
formed volume. These external forces will be, general, functions the 
spatial coordinates and time. 

consider again some deformed state the body and decompose the 
boundary traction this state into Cartesian components fn, parallel 
the fixed z-axes, respectively. The work done these forces, due in- 
crements 5ux, the displacement components, all generality 


The integral has tobe extended over the whole surface bounding the body the 
initial state. 


the case plate, the expressionfor the work the surface traction is, 
therefore, 


z=h/2 


h/2 


the preceding section, the line integral taken around the total 
undeformed cylindrical boundary the plate, while the subscripts and in- 
dicate components referred coordinates measured outwardly normal and 
along the undeformed boundary, respectively. 

Substitution the approximations for the displacements ux, 
results approximate expression for the increment the work 
done the boundary tractions, which, employing the notations 


z=h/2 z=h/2 
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h/2 h/2 
h/2 h/2 
-h/2 -h/2 


may written 


* * 
order calculate the work done body forces, let be, designate 
the components the deformed position, the body force per unit or- 
iginal volume. The work done these forces is, general, 
Substituting for the displacements ux, uy, the approximations Ux, Uy, Uz, and 


integrating through the thickness, the approximate expression for the 
work done body forces becomes 


h/2 h/2 


obtaining Eq. was assumed, for simplicity, that 


h/2 h/2 


KINETIC ENERGY 


The kinetic energy body occupying, before deformation, volume 
and having mass density given 


The dot indicates differentiation with respect time. Using for the velocities 
expressions based the approximate displacements (Eq. 3), assuming that 


4 
“4 
(11) 
: 
where 


PLATES UNDER STRESS 


independent and integrating through the thickness, the plate, ap- 
proximate expression for the increment the kinetic energy obtained 


Integrating parts with respect time and assuming that the initial veloci- 
ties the displacement components are specified, obtain 


NON-LINEAR PLATE THEORY 


Hamilton’s principle applied derive the equations motion and the 
boundary conditions. The principle states that for arbitrary time interval 
tg, the kinetic potential has satisfy 


the present case, since partial integration with respect time has already 
been carried out, 


The variation the kinetic potential must vanish for arbitrary values 
5v, and soforth. Thus, obtain from the integrand the double inte- 
grals (Eqs. 12, and 18) the following five equations motion 
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From the line integrals obtain the following boundary conditions 
du, 


PLATE UNDER INITIAL STRESS 


assumed that the state initial stress given the Cartesian com- 
ponents, which form symmetric tensor. The plate deforms its final 
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configuration, specified the displacement components ux, uy, uz, due ad- 
ditional stresses designated additional stresses are assumed 
related the displacements accordance with Hooke’s law for three- 
dimensional isotropic solid. indicated Biot (1), and also Pfliger 
(12), this approximation which introduces error the order mag- 
nitude the ratio initial stress the elastic modulus. This error due 
the fact that the deformations associated with initial stresses are neglected, 
which make the body both heterogeneous and anisotropic (12) its response 
the additional stresses. order introduce the above stresses into the non- 
linear equations motion (Eq. 21), may approximate 


o = 
may approximate 


1+E 1+E 
T =f = = 


Both relations, indicated (12), have their advantages and 
drawbacks, which briefly may summarized follows: Fig. has the 
disadvantage not describing satisfactorily the case uniaxial initial and 
additional tension since they reduce Txx Ex) instead the 
anticipated txx. Eq. 24, the tensor unsymmetrical, re- 
sulting non-symmetric equations motion. Furthermore, some terms 
the equations motion are dependent upon the arbitrary choice specifying 
Txy Sxy and forth, instead Tyx Syx and forth. 

These deficiencies can circumvented using only the relations 


Eq. 24, combined with the last three Eq. 23. Thus, taking 
obtain 


=§ +t 


Defining the initial plate stresses 


h/2 


h/2 h/2 
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and the additional plate stresses 
h/2 h/2 h/2 


and using the relations Eq. have 


noted that, depending their physical nature, the initial surface trac- 
tions and body forces may change themselves the course deformation. 
The following notation is, therefore, introduced 


x x x nn nn nn nn 


n n n n 
(29) 
x x x nn nn nn nn 
sa i i a it i i a 


which the x-component the total surface traction per unit original 
area the deformed state. the x-component the additional surface 
traction which induces the deformation. the x-component the initial 
surface traction the undeformed state. the x-component the 
change due the deformation expressed per unit original area. 
the n-component the total traction the deformed cylindrical surface per 


i 
. 
as 
ia 
; 
a 
xy xy xy Se 
4 

. 
= 


PLATES UNDER STRESS 


unit undeformed area, being the outward normal the undeformed boundary; 
the n-component the initial surface traction the undeformed state; 


the n-component the change due thedeformation, per unit 


original area; the n-component the additional traction 
cylindrical surface. 
For the body forces may put, analogously, 


x x x 


Zz Zz z 


Since initially, the plate assumed equilibrium, 


the initial stresses 
must satisfy the following equilibrium equations 


and the following boundary conditions 
nei 
nn nn 


ns ns 


tes 
+ + F +B =0 
x y 1 1 (31) 
an +x +q + B =0 
x 1 1 
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substituting Eqs. 28, 29, into Eqs. and 22, and taking into account 
Eqs. and 32, obtain the five equations motion, well the boundary 
conditions plate under initial stress. These equations will linearized 
disregarding all non-linear terms involving the additional stresses. There- 
by, are postulating that the unit elongations, shears, and rotations associ- 
ated with thedisplacements ux, Uy, Uz, are small compared unity. There- 
fore, disregard the terms involving products the derivatives the com- 
ponents the displacement and the additional stresses, for they are negligible 
comparison with the additional stresses. However, the initial stresses are 
assumed higher order magnitude additional stresses. Con- 
sequently, higher order terms involving initial stresses may the same 
order magnitude than the additional stresses. Consequently, higher order 
terms involving initial stresses may the same order magnitude the 
additional stresses and shall, therefore, retained. 
The linearized equations motion for plate under initial stress are 
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The boundary conditions are 


nn nn nn 


ns ns ns 


nn 


might noted, that identical equations motion including the boxed 
terms could obtained appropriate integration the three-dimensional 
equations motion that include initial stresses, (2). the sequel, the boxed 


terms, however, will not retained for reasons given the beginning this 
section. 


STRESS-DISPLACEMENT RELATIONS 


The relations between the additional plate stresses and the displacements 
will taken the same the absence initial stresses. mentioned 
earlier, this procedure introduces error the order the ratio in- 
itial stress component elastic modulus. 

For isotropic material the pertinent stress-displacement relations (11) 
are 
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Young’s modulus, the shear modulus, Poisson’s ratio, and fac- 
tor used for adjustment. Its value can taken equal (8) and (11). 

Relations Eq. are readily derived from Hooke’s law for isotropic 
three-dimensional solid, eliminating the z-component strain, integrating 
over the thickness coordinate and setting the integrals involving tzz equal 
zero. 


DISPLACEMENT EQUATIONS MOTION 


The stress equations motion (Eq. 33) can now converted the follow- 
ing displacement equations motion: 
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isotropic plate without initial stresses, the face-extensional motion 
(involving and v), and the flexural motion are not coupled, while aniso- 
tropic plate, such crystal, they are coupled. noted here that the in- 
itial moments and transverse shear forces couple the u-equation the 
equation, and the v-equation the equation. contrast, initial mem- 


brane forces Ni_ and Ni_ not introduce any coupling. Thus, the initial 


stresses make the plate, certain sense, anisotropic. Inasmuch the in- 
itial stresses depend, general, the space coordinates, the plate also be- 
comes heterogeneous. 


CLASSICAL PLATE THEORY THE PRESENCE INITIAL STRESS 


The passage from the non-linear plate theory discussed preceding sec- 
tions the classical non-linear plate theory (absence rotatory inertia and 
shear effects) can accomplished follows: Disregarding the rotatory in- 
ertia terms, (h3/12) and the last two Eqs. 21, the plate 
shear and may obtained terms the other plate stresses, and the 
resulting expressions may introduced into the remaining equations mo- 
tion; thus, set three equations obtained that does not involve the shear 
stresses. 
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Since the effect shear disregarded, may find setting and 


Substituting the above relations Eq. 21, set three equations motion 
obtained which excludes the effect transverse shear deformation and rota- 
tory inertia. noted that the second order terms Eq. contribute only 
higher order terms the previously-mentioned equations and can disre- 
garded. 

The same equations could have been obtained using Eq. and setting 
€yz Odirectly variation the strain 4). this case, 
however, the higher order terms Eq. cannot disregarded. 


These equations, linearized, for the case plate under initial stress, 
the manner discussed previously, are 
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The boundary conditions are 


a 
4 + — M + 


The equations motion reveal that membrane stresses not introduce any 
coupling between compressional and flexural motions, while initial 
introduce such coupling. 


Equations motion (Eq. 38) could also have been obtained disregarding 
the inertia terms and, p(h3/12) the last two Eq. 33, and 


introducing the obtained expressions for and into the remaining equa- 


spectively. 
Using Eqs. the stress displacement relations Eq. reduce 


and the displacement equations motion can readily obtained the sub- 
stitution these relations into the stress equations motion (Eq. 38). 


APPLICATIONS 


order study the influence initial stresses the frequencies the 
modes motion considered this investigation (Fig. 2), series particu- 
lar solutions the newly derived equations presented and discussed this 
section. Fig. the modes are follows: flexural, lower face- 
shear, longitudinal, upper face-shear, and thickness shear. 

Initial Normal Stress.—As first example, consider the case initial 
constant normal stress Nix well known that initial tension raises, 
while initial compression lowers the frequency the flexural vibrations. For 
the sake completeness, however, may stated that the case cylin- 


drical vibrations the flexural frequency calculated from 
the equation 


which obtained from the set Eqs. 38. 
Letting 


- 
. 
/ 
4 
| 
at, 
4 
ge 
| 
a: 


EM 3 85 
i 
ef 
3 
== 
i 
— 


June, 1960 


results 
which 
phL 


the frequency flexural vibration for small h/L. The relative effect in- 
itial compression (-N) illustrated Fig. apparent that this effect 
becomes very large for small values h/L. 


inspection Eqs. concluded that the initial normal stress does 
not effect the frequencies the longitudinal mode, and the upper face-shear 
mode. The frequency the lower face-shear mode increases due initial 
tension and decreases due initial compression; this effect is, general, 
exceedingly small. 

The influence initial normal stress the frequencies the predomi- 

nantly shear mode has been investigated (15). was found that this frequency 
increases initial tension and decreases initial compression. However 
this effect general very small. For example for 10-3 and h/L 
the frequency will change amount the order magnitude 
10-7. Notice that since the derived equations motion are valid only the 
elastic range, the value for most engineering materials could not exceed 
the order magnitude 10-3 

Initial Bending Moment.—As second example, study the influence, 
frequencies vibration, constant initial moment For face shear 


apparent that the initial bending moment introduces coupling between 
the upper and lower face-shear modes. 
seek solution the above equations the form 


v2 V sin e . . » (46a) 
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and obtain the frequency equation 


which 


Noting that, order that the stress does not exceed the elastic limit, must 


not exceed the order magnitude conclude that the frequency 
the upper face-shear mode increases due the initial moment. This effect 


increases with h/L, but is, general, exceedingly small. also observe 


FIG, 4.—PLATE SUBJECTED INITIAL TRANSVERSE SHEAR 


that the frequency the lower face-shear mode decreases due the initial 
moment. This effect also, increases with h/L but generally very small. 

Initial Shear.—As alast example, consider the influence, the frequen- 
cies free vibration, uniform initial shear force, seen 
from the equations equilibrium for the initial state stress (Eq. 31), 
state uniform initial shear cannot produced edge loading, but only 
external bending moment mi, uniformly distributed over the plate (Fig.4). 
This moment may result surface uniformly distributed 
over the two plate faces. the plate deforms, the applied moment mi, the 
surface shear tractions constituting this moment, could change both magni- 
tude and direction, depending their physical nature. 


the frequency the thickness-shear vibration. The solution Eq. 
¥ 
x = 
Q 
. 
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For this example will assumed that the magnitude the applied shear 
tractions could change during deformation such manner that for any ele- 
ment the resulting bending moment remains constant. canbe shown that this 
force system conservative. 

alternative, the magnitude the surface could remain 
constant during deformation. this case, however, the force system non- 
conservative. more specifically, Ziegler’s (13) terminology, 
tory.” 

regard the deformed direction the applied surface shear tractions, 
could assume that the resulting shear forces either rotate withthe element 
not change direction after deformation. The resulting applied bending 
moment, within the accuracy the present theory, will the same both 
cases. 


For the above mentioned conservative loading, and =0. Con- 
sidering only free cylindrical vibrations 0), the equations 
motion (Eq. 38) reduce 


The extensional motions are coupled the flexural motions or- 
der study the effect the flexural vibrations, first approxima- 
tion, the compressional inertia neglected. The equation then becomes 


the well known equation for cylindrical motions plate under initial com- 
pression Thus, may anticipated that the initial shear force will have 
pronounced effect the flexural frequencies. Furthermore, noticed 


that the inertia term neglected, Eq.51 reduced the buckling equation 
the Euler column. 


Assuming the solution Eq. the form 


obtain the frequency equation 
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the vibrations are remain finite, all the values must real. 
Therefore, buckling occurs for 


For the stresses the plate not exceed the elastic limit, must not 
larger than the order magnitude 10-3 Eh. For the sake simplicity, 
choose 0.305 and 10-3 and obtain 


From this relation apparent that the initial decreases 
the frequency the flexural mode. The relative effect becomes large for 
small value h/L; this illustrated which, from the above equa- 
tion, the ratio plotted function h/L. 

For the chosen value the plate will buckle h/L= 10-3. This indicates 
that plate must have very small h/L ratio order buckle the elastic 
range, due constant transverse shear. 


the longitudinal inertia not neglected, the frequency equation 


the frequency longitudinal vibrations the absence initial stresses, 
and given Eq. 44. 
The two roots the Eq. are 


1/2 


are small compared unity, the second term under the square root the 
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above expression small compared one, and binomial expansion results 


+ W 
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The higher frequency can then written approximately, 


seen that the initial shear stress raises the higher (longitudinal) frequency, 
but only small amount. This effect increases with h/L but is, general, 
exceedingly small. The lower (flexural) frequency 


The second term the parenthesis above represents the influence longi- 
tudinal inertia and seen small compared unity. this term 
dropped, the expression for becomes identical the one given Eq. 55. 
evidenced from this example, the initial shear stress has large effect 
only the flexural frequencies vibrationfor very slender plates(h/L <1). 
The effect initial shear the longitudinal frequency exceedingly small. 
the case face shear vibrations, the frequency equation for motion un- 
der the influence initial shear analogous the one obtained for motion 
under the influence initial bending moment. may, therefore, conclude 
that initial shear increases the frequency the upper face-shear mode, and 
decreases the lower. These effects increase withh/L but are, 
general, exceedingly small. 
The influence initial transverse shear stress the frequency the 
predominantly thickness-shear mode has been investigated (Eq. 15). was 
found that initial shear stress very small amount, 
even for relatively large values h/L. 


CONCLUSIONS 


the basis the results obtained applying the equations motion de- 
rived this paper number specific cases, may conclude that the ef- 
fect initial stresses the frequencies the two face shear modes and 
the predominantly compressional and thickness shear modes exceedingly 
small. The relative effect initial membrane stresses (Nix and 


transverse shear stress the frequencies the predominantly flexural 
mode increases h/L decreases. the elastic region, however, the effect 
initial membrane stresses becomes large range h/L that has great 
practical importance, while the effect initial transverse shear becomes 
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significant for values h/L too small great practical importance. 
For instance, the frequency the predominantly flexural mode changes 
50% for plate with h/L 0.03 when under the influence initial tension 
magnitude 10-3 Eh. produce the same change when plate subjected 
initial transverse shear the same magnitude, the ratio h/L has ap- 
proximately 0.0017. 

The qualitative influence uniform initial stresses, the frequencies 
cylindrical vibration, summarized Table 


TABLE 1.—QUALITATIVE INFLUENCE UNIFORM INITIAL STRESS 
THE FREQUENCIES CYLINDRICAL VIBRATIONS 


Initial Flexural Longitudinal Upper Thickness 
Stress face-shear shear 


Tension increases increases 


decreases decreases increases increases increases 


increases 
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APPENDIX NOTATION 


List symbols used this paper: 
Final Cartesian coordinates element. 


spectively. 


Plate thickness. 


the z-x and z-y planes, respectively. 


Total strain energy. 


(i, Trefftz’s components stress per unit original area 


(i, Components stress per unit original area (Fig. 1). 


z-direction, respectively. 


Strain energy density. 


- 
| 
— 
Moment decreases increases 
: 
Sar, 
4 
4 


Plate stresses defined Eq. 

Kxx, Kyy, Kxy 

Coordinates measured along the outward normal and 
along the tangent the undeformed cylindrical bound- 
ary respectively. 

Work done boundary tractions. 

Surface body before deformation. 


Plate area before deformation. 


Boundary plate forces and moments defined Eq. 11. 


be, Components body force. 

Work done body forces. 

Bn, Components body force defined Eq. 15. 
Kinetic energy. 

Mass density. 

Time. 

Kinetic potential. 


initial boundary plate force 
Plate compressional modulus. 
Plate flexural modulus. 


Young’s modulus. 
Poisson’s ratio. 
Shear modulus. 


Adjustment factor. 


Circular frequency the absence initial stress. 
Frequency the first thickness-shear mode the ab- 
sence initial stress. 


Frequency the flexural mode for large wave-lengths 
the absence initial stress. 


Frequency the compressional mode forlarge wave- 
lengths the absence initial stress. 
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Superscripts: 


approximate 
additional 
initial 


the boundary 
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ULTIMATE STRENGTH OVER-REINFORCED BEAMS 


SYNOPSIS 


The flexural strength rectangular reinforced concrete beams controlled 
compression investigated analytically, using stress distribution derived 
from experimental stress-strain curves. Empirical expression for the ulti- 
mate bending moment presented and the calculated moments fifty-nine 


test beams that failed compression are compared with the ultimate test mo- 
ments. 


INTRODUCTION 


American design practice, the ultimate flexural strength reinforced 
concrete beams determined means empirical design formulas based 
the results numerous For beams controlled yielding the 


Note.—Discussion open until November 1960. extend the closingdate one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Graduate Student, Civ. Engrg. Dept., Northwestern Univ., Evanston, and Assoc. 
Development Engr., Structural Development Sect., Portland Cement Assn. Labs., Skokie, 

Assoc. Prof. Civ. Engrg., Northwestern Univ., Evanston, 

“Ultimate Strength Reinforced Concrete American Design Practice,” 
Hognestad, Proceedings, Symposium the Strength Concrete Structures, London, 
May, 1956. Also, Portland Cement Assn. Development Dept. Bulletin D12. 


“Report ASCE-ACI Joint Committee Ultimate Strength Design,” Proceedings 
ASCE, Vol. 81, October, 1955. 
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tensile reinforcement these formulas, which were also derived 
give close approximation the ultimate strength. The flexural strength 
beams controlled crushing concrete may estimated, assuming that the 
strength the beams proportional tothe concrete cylinder strength.4,6 Be- 
cause the undesirability compression failures, formulas based this 
assumption are designed give conservative value for the calculated strength. 

When greater accuracy the calculated strength beams controlled 
compression required, other methcds may used. analysis based 
assumed rectangular stress distribution and ultimate concrete strain 
gives satisfactory the present paper the flexural strength rec- 
tangular beams controlled compression investigated analytically, using 
stress distribution derived from experimental stress-strain curves. 


CONCRETE STRESS DISTRIBUTION FLEXURE 


The stress-strain relationship concrete flexure was determined ex- 
perimentally for concrete various strength the ages 14, 28, and 
For the purpose this investigation, these stress-strain curves were 
approximated second-degree curves having the general equation 


which denotes the stress, the strain, and and are coefficients 
which are functions Solving for the stress yields 


The coefficients Eq. solving four simultaneous equations, 
obtained substituting Eq. four given values and the corre- 
sponding values stress The values strain used, inches per inch, were 
0.0006, 0.0012, 0.0018, and 0.0024. 

The values stress corresponding each strain were first ex- 
pressed function the cylinder strength means the equation 


The coefficients and Eq. which are functions were deter- 
mined from the twenty stress-strain curves given Hognestad, al,8 
the method least squares, and are presented Table 

With the coefficients given Table the values stress corresponding 
thefour given values strain were determined for concrete strength 
ing 1,000-psi increments 6,000 psi. Then the coefficients Eq. 


“Ultimate Strength Criteria for Reinforced Concrete,” Kriz, Proceedings, 
Vol. 85, No. July, 1959. 
“Plastic Theory Reinforced Concrete Design,” Whitney, Transactions, 
ASCE, Vol. 107, 1942, pp. 251-282. 
“Rectangular Concrete Stress Distribution Ultimate Strength Design,” 
Mattock, Kriz, Hognestad and “Rectangular Concrete Stress Distribution 
Ultimate Strength Design,” published. 
“Concrete Stress Distribution Ultimate Strength Design,” Hognestad, 
Hanson and McHenry, Proceedings, ACI, Vol. 52, 1955, pp. 455-479. Also, Portland 
Cement Assn. Development Dept. Bulletin D6. 
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were calculated, and the results are given Table The strains are 
given 10-3 inches per inch, and the stresses kips per square inch. 
The stress-strain curves obtained this manner are shown Fig. 


TABLE 


0.1134 
0.04016 -0.04800 


0.05359 -0.01559 
-0.02019 


TABLE 


-0.07042 


The following properties the stress-strain curves were determined for 


later use numerical integration using the and strain incre- 
ments 0.0002 in. per in.: 


ULTIMATE STRENGTH ANALYSIS RECTANGULAR BEAMS 


The ultimate flexural strength rectangular reinforced concrete beams, 
shown Fig. was determined computing the resisting moment the 
beams, assuming increasing values the strain the extreme compressive 


(1) (3) (3) (4) 

0.0012 

0.0024 

(2) (3) (4) (5) 
1,000 -1.262 1.746 
3,000 -3.434 -6.751 22.58 
e 

4,000 -6.603 -14.26 56.99 

5,000 -1.196 -31.41 137.8 

and 
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fibers with 0.004 in. per in. maximum. The highest value the re- 
sisting moment thus obtained was designated My, the ultimate bending moment. 
The following assumptions were made the computations: 


(a) Plane sections normal the axis remain plane after bending. 

(b) Tensile strength concrete may neglected. 

(d) The modulus elasticity the reinforcement 106 psi. 


(e) The stress the reinforcement constant after the yield stress 
reached. 


The yield strengthof the reinforcement was assumed 60,000 psi, the 


maximum stress recommended the ACI Building Code (ACI 318-56), Section 
A603. 


FIG, 2.—CONDITIONS REINFORCED CONCRETE BEAM BENDING 


With these assumptions, the equations equilibrium and compatibility 
strains give the resisting moment 


which 
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and 


means through the properties the approximate stress- 
strain curves, the ultimate strength flexure was determined for beams with 
percentage tensile reinforcement varying increments 0.005 from 
0.005 0.050, and concrete compressive strength varying incre- 
ments 1,000 psi from 1,000 psi 6,000 psi. The assumed strain 
the extreme compressive fibers was varied increments 0.0002 in. per 
in. until the maximum moment strain 0.004 in. per in. was reached. 
All computations outlined previously were programmed for LGP-30 elec- 
tronic computer. 

concrete compressive strengt with the percentage reinforcement 
parameter. The area the graph divided heavy line into two regions: 
(a) one below the heavy line, which the dashed lines represent the ultimate 
bending strength beams controlled tension, and (b) the one above the di- 
vtding line, which the full lines represent the strength beams controlled 
compression. The dividing line was determined the intersections the 
two systems lines, and, thus, represents the ultimate strength beams with 
balanced reinforcement for the yield stress 60,000 psi. The ultimate 
strength beams controlled tension, obtained the present method, 
close agreement with the values obtained Eq. the ACI Building Code 
(ACI 318-56). 

The variation the ratio My/ suggests the possibility approxi- 
mating the ultimate flexural strength rectangular reinforced concrete beams 
controlled compression means the expression 


or, specifically, 


Eq. 11b was used compute the ultimate strength fifty-nine test beams 
which failed The results are presented Table 
and Fig. The distribution the ratio between the test and the calculated 


“Compressive Strength Concrete Flexure Determined from Tests Rein- 
forced Beams,” Slater and Lyse, Proceedings, ACI, Vol. 26, 1930, pp. 831- 
874, 

“Comparative Tests Concrete Beams Reinforced with Isteg and Hot Rolled De- 
formed Bars,” Columbia University, Report No. 2507, 1941. 

Reinforced Concrete Beams with Recommendations for Attaining Balanced 
Design,” Cox, Proceedings, ACI, Vol. 38, 1942, pp. 65-80. 

“The Ultimate Strength Reinforced Concrete Beams,” Lash and 

Brison, Proceedings, ACI, Vol. 46, 1950, pp. 457-470. 
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TABLE 3.—COMPRESSION FAILURES REINFORCED CONCRETE BEAMS 
REINFORCED TENSION 


1 
Slater 
and 
Lyse 


Columbia 
University 


w 


oom eo 


wo 


Cox 


Brison 


Average, all beams 
Standard Deviation, all beams 


+ 


ber Test Calculated 
(1) (2) (3) (4) (5) (6) (7) 
2,820 64.8 0.030 0.346 0.372 
10A 3,810 64.8 0.040 0.354 0.344 
3,550 61.37 0.0341 0.406 0.346 
Cll 3,550 62.60 0.0345 0.394 0.347 
3.550 63.04 0.0334 0.386 0.345 
C12 3,550 0.0328 0.365 0.344 
122 1,700 53.4 0.0176 405 0.402 
123 1,700 53.4 0.0264 0.417 
124 1,700 53.4 0.0352 0.431 
125 1,700 53.4 0.0440 0.446 
142 1,700 0.0244 0.414 
143 1,700 48.1 0.0368 0.434 
144 1,700 48.1 0.0488 0.454 
224 3,100 53.4 0.0352 0.368 
225 3,100 53.4 0.0440 0.383 
243 3,100 48.1 0.0368 0.371 
244 3,100 48.1 0.0488 0.390 
214 3,100 55.2 0.0308 0.374 0.361 
215 3,100 55.2 0.0388 0.372 0.374 
235 3,100 48.1 0.0388 0.381 0.374 
253 3,100 50.6 0.0520 0.388 0.395 
325 4,500 53.4 0.0440 0.320 0.320 
343 4,500 48.1 0.0368 0.297 0.308 
344 4,500 0.0488 0.342 0.328 
425 5,800 53.4 0.0440 0.293 0.261 
444 5,800 48.1 0.0488 0.300 0.269 
4205 1,970 39.2 0.0400 0.449 0.427 
4206 1,930 44.4 0.0475 0.464 0.441 
4308 3,330 42.8 0.0454 0.367 0.375 
4407 4,170 40.8 0.0367 0.305 0.322 
4408 4,490 43.4 0.0471 0.346 0.325 
6203 2,150 88.0 0.0140 0.346 0.376 
6204 2,150 75.8 0.0200 0.333 0.386 
6205 1,950 75.8 0.0225 0.415 0.399 
6207 1,915 75.2 0.0385 0.391 0.427 
6208 2,120 75.2 0.0391 0.407 0.419 
6303 3,290 72.0 0.0147 0.270 0.326 
6304 2,760 75.8 0.0233 0.347 0.364 
6305 3,200 74.0 0.0286 0.338 0.352 
6306 2,760 75.2 0.0394 0.359 0.390 
6404 4,490 75.8 0.0226 0.263 0.285 
6405 4,140 74.0 0.0280 0.293 0.310 
6406 4,190 75.2 0.0390 0.317 0.325 
6407 4,190 62.1 0.0408 0.314 0.328 
6504 4,870 75.8 0.0233 0.247 0.269 
6505 4,450 65.0 0.0371 0.299 0.310 
6506 5,450 75.8 0.0458 0.292 0.280 1,0 


OVER-REINFORCED BEAMS 


1000 2000 3000 4000 5000 6000 
Cylinder (psi) 


FIFTY NINE BEAMS CONTROLLED COMPRESSION 


moment shown histogram, Fig. The average value this ratio for 
the fifty-nine beams 1.007, with standard deviation 0.076. 


CONCLUSIONS 


The ultimate strength rectangular reinforced concrete beams, when con- 
trolled compression, depends primarily concrete strength and, less- 
degree, the percentage tensile reinforcement. The relationship be- 
tween the ultimate strengthin bending and the two factors nonlinear. However, 
the ratio may assumed decrease linearly with increasing 
concrete strength (within the practical range 2,000 psi 6,000 psi) for con- 
stant percentage reinforcement. Similarly, the ratio given 
concrete may assumed proportional the percentage re- 
inforcemént. equation basedon these two assumptions gives close approxi- 
mation of, the ultimate bending moment, confirmed comparison test 
values with the computed values ultimate moment fifty-nine test beams, 
which faijed compression. 
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FIG, 5.—HISTOGRAM FOR COMPRESSION FAILURES SINGLY 
REINFORCED REINFORCED CONCRETE BEAMS 
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NOTATION 
width rectangular section; 
distance extreme compressive fiber centroid tensile re- 
inforcement; 
modulus elasticity tensile reinforcement; 
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concrete cylinder strength time test; 
stress tensile reinforcement; 
yield point reinforcement; 


ratio distance between extreme compressive fiber and neutral 
axis distance between extreme fiber and centroid tensile 


ratio distance between extreme compressive fiber and resul- 
tant compressive stresses distance between extreme fiber 
and neutral axis, defined Eq. 10; 


resisting moment; 


ultimate moment; 


percentage tensile reinforcement defined which 
denotes area tensile reinforcement; 


constant coefficients; 
strain concrete; 
strain extreme compressive fiber; and 


strain tensile reinforcement. 
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EXPERIMENTAL STUDY BEAMS ELASTIC FOUNDATIONS 


Robert Thoms! 


SYNOPSIS 


The design and use structural model for the solution shear and mo- 
ment beam elastic foundation presented. Through use the 


model, typical problem laterally loaded pile foundation solved 
example. 


THE PROBLEM 


Shear and moment loaded beam elastic foundation can calcu- 
lated the solution the differential equation relating beam deflection and 
load available and the deflection-reaction characteristics the support- 
ing foundation are known. will assumed this study that deflection- 
reaction characteristics the supporting foundation are known. The differ- 
ential equation relating load and the form, 


(1) 
Numerical methods solving Eq. for the bending beam continuous 


elastic foundation generally make use anassumption Winkler, that the re- 
action forces the foundation any point are proportional the deflection 
the beam that same The effect this assumption represent 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 
Inst. Engrg. Mechanics, The Univ. Texas, Austin, Texas. 
“Strength Materials, Part II,” Timoshenko, New York, Van Nostrand, 


107 


Sere 
= 
Le 
‘ 
iy. 
: 


108 June, 1960 


the beam being supported closely spaced coil springs acting independent- 
each other. The solution foundation problem based the Winkler 
assumption neglects any existing shear between adjacent supports. also 
assumed that negative reactions may exist the elastic beam. 

The method solution presented herein alternate procedure nu- 
merical methods for finding deflections beam elastic foundation. 
This method has the advantage being more rapid than numerical computa- 
tion means desk computers, and less expensive than computation involv- 
ing large electronic computers. Foundation problems with any pattern vari- 
ation reaction modulus with distance along the beam can solved. 


SMOKING 


FIG, 1.—MODEL LOADED 


more, this method can readily accommodate changes the elastic beam, 
and could extended provide solutions where structural restraints are in- 
volved the beam. Although this method has the property being less exact 
than numerical methods, should noted than the accuracy with which soil 


reaction characteristics can determined does not necessitate more exact 
solution. 


METHOD SOLUTION 


The method used involves structural model that physical adaptation 
the Winkler beam. Basically, the model consists anelastic beam support- 
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closely spaced coil springs. The model was designed for the purpose 
predicting the deflection curve loaded elastic beam elastic foun- 
dation. Reactions the beam may found use the predicted deflection 
curve the beam and the reaction modulus the supporting foundation. Shear 
and moment the beam may then calculated successive integration 
the reactions the beam. 

problem with supporting foundation possessing non-linear deflection- 
reaction characteristics can solved use iteration process. The 
iteration process necessary for the model, essentially linear system, 
represent non-linear effects. 

the example problem this study the model was designed for the purpose 
predicting the deflection curve laterally loaded pile The 
predicted deflections the pile, obtained use the model, were used 
iteration process effective secant soil reaction moduli, 
the supporting soil. The iteration process was used adjusting the model 
represent the non-linear effects the supporting soil reaction the pile. 
Supporting soil reactions the pile were then obtained multiplying the pre- 
dicted pile deflections the determined effective secant soil reaction moduli. 
With the reactions the pile available, complete solution the problem 
was possible. 

The method solution described for the example pile problem different 
from existing numerical methods only that the model was used obtain the 
deflection curve the loaded pile. The advantage the model that may 
used obtain the deflection curve quite rapidly, thereby greatly reducing 
the time spent the iteration process used determining the effective secant 
soil reaction moduli, 


SIMILITUDE REQUIREMENTS 


The system represented the model consists (1) the elastic beam 
and (2) the supporting foundation. The elastic beam was represented the 
model cylindrical steel rod with flexural rigidity and length determined 
flexural similitude requirements. The supporting foundation was represent- 
the model series closely spaced elastic supports. Each support 
viding there are model supports equally spaced, each model support repre- 
sents increment supporting foundation length, 

required that the model have deflection curve similar that the 
actual pile prototype undergoing study. this, the model must 
satisfy certain similitude requirements. These requirements may deter- 
mined the methods dimensional analysis. 

The deflection distance from the loaded end the beam, Fig. 


(Ey and (L)p, the prototype factors. 


analysis: 


Engineering,” Glen Murphy, New York, The Ronald Press Co., 
1950, pp. 112-113, 118-123. 
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and 


required that Qm, therefore 


and 


Eqs. and may combined give, 


Eq. was the basic relation observed, along with Eq. design and ad- 
justment the model. should remarked that the structural model referred 
this study not generally simple scaled-down version prototype 
beam, this would have limited its versatility. The basic requirements, Eqs. 


and were combined into Eq. permit convenient representation any 
prototype the model. 


PREDICTION EQUATIONS 


The relation between deflection the model and prototype may estab- 
lishedon the basis dimensional analysis. This relationis generally referred 
the prediction equation for the 

Since 


and Eq. must true, the combination Eqs. and yields the prediction 
equation for the model: 


Ibid., 99. 
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similar fashion, the prediction equation for the model loaded with 


applied concentrated moment is, 

Eq. may written 


DESCRIPTION THE MODEL 


The model used for this study consists cylindrical steel rod 
ing beam and series adjustable coil spring supports representing the 


L 
FIG, 2.—TYPICAL BEAM 


foundation support. Fig. shows the model with loadapplied one end. The 
elastic rod seen the center the supporting frame the model represents 
the beam. This particular rod in. length. Outside dimensions the 
relatively rigid supporting frame are in. Each coil spring con- 
nected the supporting frame the model with device means which 
the length the coil spring supporting the elastic rod can quickly and con- 
veniently changed, thereby allowing any support range spring constants. 

The model may loaded weight suspended over pulleyor free- 
hanging weight. Care must that closureof any coil springs avoid- 
ed. Deflections the rod are obtained recording, both before and after 
loading, the vertical positions pointers clamped the rod. eliminate 

reading parallax, engineers transit was used approximately from 
the model note the positions the pointers the grid immediately behind 
the rod. 


1 
a 
4 


112 June, 1960 
USE MODEL PROBLEM SOLUTIONS 


The model may used solve for the deflection beam elastic 
foundation possessing (a) linear deflection-reaction characteristics, (b) non- 
linear deflection-reaction characteristics. 

The following procedure used solve the linear case: 


(1) Preliminary calculations are used selecting suitable elastic rod 
represent the beam. Prototype characteristics and initial arbitrarily fixed 
model characteristics are substituted Eq. determine the approximate 
dimensions satisfactory rod. 

(2) the basis the prototype conditions and the selected elastic rod, the 
coil springs may adjusted represent the foundation reaction moduli. 

(3) The model loaded with any convenient magnitude load that does not 
cause closure the coil springs, and deflections measured. 

(4) The deflections may used predict prototype deflections through use 
prediction equations. 

(5) The predicted deflections the prototype may used find the foun- 
dation reactions the elastic beams through use the reaction moduli. 

(6) The moment and shear diagrams can obtained successive inte- 


gration the reactions the beam, and then used calculate stress the 
beam. 


The non-linear case solved with the model the following procedure: 


(a) This step identical step for the case. 

(b) the basis the prototype conditions and the selected elastic rod, 
estimates deflection can made and corresponding secant moduli obtained 
from the available foundation reaction curves. The coil springs may then 
adjusted represent these initial estimates reaction moduli. 

(c) The model loaded convenient magnitude load that does not 
cause closure the coil springs, and deflections measured. 

(d) The deflections may used predict prototype deflections through use 
prediction equations. The predicted prototype deflections are plotted the 
foundation reaction curves and new estimates the reaction moduli determined. 

(e) The model supports are adjusted torepresent the secondestimate the 
effective foundation reaction moduli. 

Steps and are repeated until the effective secant foundation reac- 
tion moduli are determined completely. The effective secant foundation reac- 
tion moduli are determined when the reaction moduli obtained the course 
the trial and adjustment process converge sufficiently constant values. 

(g) The final the model with supports correctly adjusted used 
predict deflections the prototype. The predicted deflection the proto- 
type transformed into foundation reaction the pile use the effective 
secant foundation reaction moduli determined previously. 

(h) This step identical step the solution the linear case. 


EXAMPLE PROBLEM 


The following example illustrates the solution hypothetical problem 
through use the model. The hypothetical problem the analysis moment 
and shear laterally loaded pile foundation. 
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The analysis entailed designing laterally loaded pile involves two major 
problems: (1) Soil reaction curves must obtained, and (2) the deflection 
the pile must calculated order solve for shear and This 


study concerned primarily with the second problem, and soil characteristics 
are assumed known. 


There are, generally, two methods solving for the shear and moment 
laterally loaded pile: (1) methods involving limit analysis, and (2) methods 
based elastic Although physical model may adjusted 


KIPS. 


INITIAL 


EI= 


(inches) Initial 


EXAMPLE PILE HYPOTHETICAL SOIL 
REACTION CURVES 


P=50 KIPS APPLIED LATERALLY 
INCHES. 
INITIAL ASSUMED VALUES 
PER INCH FOR 200", 
PER INCH FOR 400", 
PER INCH FOR 600". 
FIG, 3.—DATA ASSUMED FOR PROTOTYPE PILE 


Solutions for Laterally-Loaded Piles With Soil Modulus As- 
sumed Proportional Depth,” Lymon Reese and Hudson Matlock, Proceedings 


Eighth Texas Conference Soil Mechanics and Foundation Engineering, 1956, pp. 
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use either approach, this study will restricted ananalysis based 
representation the pile anelastic beam anelastic foun- 
dation. 

Fig. shows the given prototype data and the hypothetical available soil re- 
action curves various depths the soil. The soil was divided into three 
layers, each layer possessing the soil reaction curve shown. the soil re- 
action curves are non-linear, follows that the procedure previously described 
for use the model representing non-linear foundation systems will used 
solving the problem. may pointed out that fifteen layers soil, 
each possessing individual soil reaction curve, could have been represented 
the model shown Fig. 

initially four variables the model that may adjusted 
senting the actual pile. These variables are: (1) the length the elastic rod 
representing the pile, (2) the spacing the supports, (3) the flexural rigidity, 
the model pile, and (4) the spring constants the supports. 


2nd. TRIAL 


ist TRIAL 


INITIAL 
=°500'bs/in 


SOIL REACTION 
CURVE 


SOIL 


FIG, 4.—SOIL REACTION CURVE 


The length the elastic rod representing the pile was arbitrarily selected 
in., thereby establishing scale ratio and model. 
The supports were arbitrarily spaced 4-in. equal intervals, with the first 
and last support being in. from the ends the elastic rod. This gave total 
fifteen supports. The length soil support the prototype represented 
each model support was then in. The remaining two variables, the flexural 
rigidity the elastic rodand the spring constants the model supports, were 
determined from Eq. with prototype data substituted, and consideration 
properties available coil springs and rods. Coil springs with range 
spring constants from approximately 0.55 per in. 7.00 per in. depend- 
ing the length coil spring active supporting load were used for the 
model supports. steel cylindrical rod, in. length, with change sec- 
1,720 for 20in. $60 in., was selected represent the proto- 
type pile. 
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With the establishment the dimensions the steel rod and the spacing 
supports, was possible use Eq. determine the spring constants 
the model supports such that the initial estimated secant soil reaction moduli 
shown Fig. would represented. After all spring constants had been de- 
termined, the model supports were adjusted possess the necessary constants 
the basis the assumption that the spring constant any coil spring was 
inversely proportional the length spring active supporting load. 

The model was loaded with force the end 
the rod corresponding the top the prototype pile. Deflections were mea- 
suredand used predict prototype deflections with Eq. 10. use the pre- 
dicted deflection the prototype and the soil reaction curves for each soil 
support, new soil reaction modulus was determined. This modulus was taken 
the slope straight iine from the origin the soil reaction curve the 
point the curve corresponding the predicted deflection. This procedure 
generally gave better estimate the effective secant soil reaction modulus 


TABLE 1.—DEFLECTIONS MODEL AND PROTOTYPE 


Location Model Measured Deflection |Location Predicted Deflection 


(2) (3) (4) 


than the initial assumed value. The model supports were then readjusted 
reproduce the effect the newer soil reaction moduli. The model was re- 
loaded and deflections recorded. Prototype deflections were predicted the 
basis the second set model deflections, and soil reaction moduliagain de- 
termined. This procedure was repeated until the secant soil reaction moduli 
converged effectively constant values. 

The soil reaction curve the soil support in. shown Fig. 
The initial and subsequent soil reaction moduli used the calculations are 
also shown. Only three trials were necessary determine the effective se- 
cant soil reaction modulus, illustrating the rapid convergence the method 
this study. 

Table shows measured deflections the model pile. Soil reactions the 
actual pile obtained soilreaction curve and the final values the pre- 
dicted deflections are shown Col. Table 


0.03 240 0.07 

0.03 280 -0.07 

0.03 320 0.06 

0.00 440 0.00 

-0.01 560 -0.01 
-0.01 600 -0.02 
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The shear and moment diagrams computed from the soil reactions for the 
laterally-loaded pile are shown Figs. and The diagrams were obtained 
integrating successively the reactions from the loaded end the pile. 

Apparent error this solutionis indicated the failure the moment and 
shear diagrams close zero 580 in. The magnitude error shown 
typical the method solution this study. may pointed out that 
the reactions were integrated successively from the loaded end the beam, 
any smallerrors shear the loaded end would have become magnified error 
moment the unloaded end. neighborhood the more important 
maximum moment, the value moment obtained the integration process 
not easily affected the moment near the unloaded end the pile. 


TABLE 2.—RESULTS FOR THE PROTOTYPE PILE 


Location Predicted Average Effective Soil Reactions 


Prototype Deflections Over Secant Soil Pile 
Soil Supports Reaction Moduli (L) 


100 
140 
180 
220 
260 
300 
340 
380 
420 


cocoooorn 


With the shear and moment the pile obtained reasonable degree 
accuracy, the stress could calculated the conventional manner. 

After some experience inuse the model, doubt closer initial estimates 
the effective secant soil reaction moduli could made, thereby reducing 
the time required obtain satisfactory solution. 


CONCLUSIONS AND FURTHER USES THE MODEL 


The major advantage the model that makes possible reasonably 
accurate solution complicated problem with small amount labor and 
expense. The model can used solve for shear and moment elastic 
beam elastic foundation with any the following characteristics: (1) any 
variation with (2) changes the beam, and (3) structural re- 
straints the beam. Shear and moment elastic beam with the foregoing 
characteristics can solved numerical computer methods. However, 
these methods are not generally simple and convenient use. 


-1.30 382 
-0.91 538 
-0.29 800 
0.07 1750 
0.06 2000 
0.04 1820 
0.02 1540 
0.01 640 
460 0.00 910 
500 0.00 850 
540 0.00 1250 
580 -0.01 1250 
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difficult problem that may solved use the model laterally 
loaded pile foundation with all the previously listed characteristics. This 
problem cannot readily analyzed non-dimensional solutions such are 
available for pile foundations with constant and 


15.2 


LAST SUPPORT 


UNITS KIPS SHEAR 


FIG, 5.—SHEAR DIAGRAM FOR TYPICAL PROBLEM 


LAST SUPPORT 


2000 


2500 


UNITS INCH-KIPS MOMENT 


—3000 


FIG, 6.—MOMENT DIAGRAM FOR TYPICAL PROBLEM 


would also possible for the model extended solve problems 
grid-type beams. this case, the model would a.more exact representa- 


tion the prototype than the case prototype involving continuous sup- 
porting foundation. 


4 
12.8 10.3 
5.5 
-10.6 
-30,2 
-50 
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benefit obtained through experience with the model would insight into 
the action the beam-foundation system. The model could also used 
supplement numerical methods. Gross errors computation could de- 


tected rapidly from rough check obtained with the model. 
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NOTATION 


reaction modulus supporting foundation, pounds per inch de- 


flection per inch beam; 


effective reaction modulus supporting foundation, pounds per 


inch deflection per inch pile; 


secant soil reaction modulus, pounds per inch deflection per 


inch pile; 


flection per inch pile; 


flexural rigidity, pound-inches2; 


length beam, inches; 

moment, inch-pounds inch-kips; 
number elastic supports; 

applied lateral load, pounds kips; 
distance along beam, inches; 

deflection, inches; 


quantity parenthesis characteristic model; and 


meee 


quantity parenthesis characteristic prototype. 


effective secant soil reaction modulus, pounds per inch de- 


spring constant elastic support, pounds per deflection; 
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ASYMMETRICAL BENDING CONICAL SHELLS 


SYNOPSIS 


The equations Goldenveizer’s theory are specialized the case con- 
ical shells. method solution developed applicable truncated cones 
finite length simply supported along two generators and having any type 


support the ends. Numerical results are obtained using electronic 
computer. 


INTRODUCTION 


Recent the problem asymmetrical bending conical 
shells have been primarily concerned with the derivation equations, while 
appears that relatively little progress has been made the direction 
finding solutions from which numerical results can readily obtained. The 
purpose this paper present practical method which such solutions 
can found, and illustrate the method with numerical examples. 

The method based well-known separation variables technique and 
the development power series form solution relative ordinary 
part the system differential equations. The series converge very slowly, 
but this difficulty overcome using high speed digital computer for their 


Note.—Discussion open until November 1960. the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Asst. Prof., Engrg. Mechanics, Univ. Kansas, Lawrence, Kans. 

“Thin Circular Conical Shells Under Arbitrary Loads,” Hoff, Journal 

Applied Mechanics, Vol. 22, No. 1955, pp. 557-562. 

Donnell Type Theory for Asymmetrical Bending and Buckling Thin Conical 
Shells,” Seide, Journal Applied Mechanics, Vol. 24, No. 1957, pp. 547-552. 
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evaluation. The method applicable cases where cross-sections are com- 
plete circles, and the special case axial symmetry, although the present 
work shall not concerned with the latter. 

Before particular problem can attacked, general theory shells 
must accepted. The question what constitutes suitable theory has been 
the source considerable disagreement the literature, but for present pur- 
poses, the point view expounded the Russian theoretician Golden- 
veizer4 shall accepted. Although has some rather appealing qualities, 
this point view seems not have attracted much attention among Western 
theoreticians (perhaps because unfortunate language barrier). pro- 
posedto apply the simplest possible theory, based the Kirchhoff hypothesis, 
for which unique theorem can proved the same manner was done 
Kirchhoff for the classical theory elasticity, and for whicha (Rayleigh- Betti) 
theorem reciprocity work can proved. The formal requirements 
such theory were first established important theoretical paper.5 

One the principal points departure inthe development Kirchhoffean 
theory shells the choice stress-strain relations, that is, relations be- 
tweenthe stress resultants andcomponents deformation. Goldenveizer con- 
tends that the most suitable choice the simplest one which satisfies the re- 
quirements ofthe previously and therefore free for- 
mal contradiction. variety stress-strain relations will fulfill these re- 
quirements, and suggested that the variant most convenient for the solu- 
tion the particular problem hand When the stress-strain re- 
lations have been chosen, basic system differential equations governing 
the displacements the middle surface shell can easily obtained from 
equations governing the equilibrium stress resultants. 


BASIC DIFFERENTIAL EQUATIONS 
Let and coordinates the middle surface shell along lines 
principal curvature having radii and the coordinate out- 


ward normal the middle surface, then the length line element space 
can expressed the form 


The stress resultants may defined terms stress follows: 


“Theory Thin Elastic Shells,” Goldenveizer, Gostekhteoretizdat, 1953. 


(in 
“On the Applicability General Theorems the Theory Elasticity Thin 
Shells,” Goldenveizer, Prit. Mat, Mekh., Vol. No. 1944, pp. 1-14 (in Rus- 

sian with English summary). 
«Theory Thin Elastic Shells,” pp. 70-71, 80-82. 
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which the thickness the shell. The positive sense the stresses 
etc., that customarily taken the theory elasticity. 

With these definitions, the equations governing the equilibrium the stress 
resultants are 


and 
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and 


Noy + Ri Ry Mo; 0 & & 4.4 (3f) 


which Xj, and are components the lateral load per unit middle 
surface area the positive x2, and directions respectively. comma 
followed subscript denotes partial differentiation with respect the in- 
dicated variable. Eq. the “sixth” equilibrium equation which expresses 
the equilibrium moments about axis normal the middle surface. 

The development thin shells based the Kirchhoff hy- 


pothesis, which equivalent assuming that the components the strain 
tensor vanish, that is, 


follows from this, without further approximation, that the remaining com- 
ponents the strain tensor can expressed the form 


and 
which 
2,1 
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with and being the displacements the middle surface the positive 
and directions respectively. 

further assumed that can neglected Hooke’s law, which then 
takes the form 


where Young’s modulus and Poisson’s ratio. 

The stress-strain relations for theory shells can now obtained 
substituting Eqs. and into Eq. and performing the indicated integrations. 
the integrations are carried out with all terms containing the factor 
ing retained, the result the Fligge-Vlasov the other hand, 
the terms z/R are neglected comparison unity, and then the integra- 
tions are carried out, the result corresponds Love’s first 


“Statik und Dynamik der Schalen,” Fligge, Julius Springer, 1934. 

Differential Equations General Theory Elastic Shells,” Vlasov, 
NACA 1241, February, 1951. 

Theory Shells,” Vlasov, Gostekhizdat, 1949 (in Russian). 

“Treatise the Mathematical Theory Elasticity,” Love, Dover, 4th 
Edition, 1927, Article 329. 
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When considered from the point view what calls “authen- 
tic” accuracy, concluded that the latter procedure should adopted. Be- 
cause the inexactness the fundamental hypothesis, one cannot expect 
improve the accuracy Kirchhoffean theory shells merely applying 
the more complicated type stress-strain relations. 

Having agreed that the approximations involved Love’s stress-strain re- 
lations and the other approximations the theory are commensurate, what 
needed now convenient modification Love’s relations eliminate the 
contradictions they imply. Applying the following system, suggested Golden- 
which satisfies the requirements internally consistent theory, 


Noi (cz + 3Ry « (8d) 

where 


These particular stress-strain relations were first written down 
paper indicates that may have important bearing the present work with 
conical shells, but spite considerable effort, the author has thus far 
been unable locate copy this country.) can seen that these rela- 
tions satisfy the sixth equilibrium Eq. identically. This necessary for the 
proof the theorems previously mentioned. 

The principal difference between the system (Eq. and Love’s first ap- 
proximation the term the expressions for and This term 
has important effecton the system differential equations solved for 
and but does not imply that there will significant difference 
the calculated values and 


“Theory Thin Elastic Shells,” 71. 

81. 

New Method Calculating Thin Shells,” Novozhilov, Akad. Nauk SSSR, 
Otd. Tekh, Nauk, Izvestiia, No. 1946, pp. 35-48 (in Russian). 

“Bending Torsion Conical Shells,” Balabukh, Tsentral’ni Aero- 
Gidrodinam. Inst., Trudy, No. 577, Moscow, 1946 (in Russian). 
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When and are eliminated from the first three equilibrium equations 
(Eq. 3), and then Eqs. and are substituted, system three differential 
equations governing and obtained. Instead writing down this sys- 
tem equations for the general case, conical shells will specialized, which 


FIG, 


the case most concerned with here. This can setting (Fig. 


=tan these substitutions, the following system basic differential 
equations for conical shells obtained from Eqs. and 
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SOLUTION BASIC DIFFERENTIAL EQUATIONS 


The problem now find solutions the system differential Eq. 10. 
Suppose that the surface loads can expressed the form 


and 
x 
s B 


and the angle indicated Fig. Solutions Eq. can now written 
the form 


and 
W(x) sin ..... 


provided and satisfy the following system ordinary differential 
equations: 
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3t2 


P 2(2 - v) h2 x2 w"' + 2(1 - 2v) h2 x w' 


which primes denote differentiation with respect 

not difficult see that solutions the form (Eq. 12) identically satis- 
the boundary conditions for simple support the generators and 
(Fig. 1). setting this form can also applied complete cones. 

the solution Eq. which satisfies some prescribed bound- 
ary conditions the ends and then the problem the bending 
truncated conical shell which simply supported along two generators 
and which under the action the loads (Eq. 11) will completely solved. 
Hence, the next task develop solutions the system ordinary differ- 
ential equations (Eq. 13). 
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while Hoff2 proposes apply the method Both these 
methods are discussed who asserts that neither useful for 
the purpose obtaining numerical results. 

The Frobenius method will briefly considered. Solutions the homo- 


geneous equations (Eq. 13) are sought relative the singular point 
the form 


oo 
and 


This will represent solution root indicial equation. The point 
said regular singular point since the degree the indicial 
equation equal the order the system differential equations, namely, 
eight. Certain the roots will differ integer, and hence the general so- 
lution Eq. logarithmic terms. order estimate the inter- 
val convergence the solutions (Eq. 14) the system (Eq. 13) 
into system eight first order differential equations setting 
homogeneous system (Eq. 13) can written the form 


The the circle convergence the solutions Eq. and hence 
Eq. determined the singularities the coefficients Certain 
these have the denominator the factor [x2 which arises from 
the coefficient the second member Eq. 13. thought be- 
ing complex variable, and the view extended that the Fuchsian 
then can concluded that convergence the series solutions 
(Eq. assured only <(2h)/(t V3). This can hardly considered 
useful range for However, not unreasonable argue that, the coef- 
ficients quantities the order (h/t)2 can neglected comparison 
this event the solution will converge for all finite But these quan- 
tities are neglected, then the differential equations will longer conform 
the requirements aninternally consistent theory. This can serious under 


«Ordinary Differential Equations,” Ince, Dover, 1956, 396. 


Pohle, Journal Applied Mechanics, Vol. 23, No. 1956, pp. 322-323. 
Ince, Chapter XV. 
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some circumstances, see, for where shown that improper 
modification the stress-strain relations frame work the Kirch- 
hoff hypothesis can lead arbitrarily large error the calculation 
long cylindrical shells. But such circumstances for conical shells need not 
examined here since the indicated quantities will not neglected the 
solutions follow. 

Instead attempting apply the method Frobenius, observe that since 
the theory valid for thin shells only, follows that solutions the differ- 
ential equations (Eq. 10) are not applicable the neighborhood the apex 
conical shell. Consequently, seems logical seek solutions the sys- 
tem (Eq. 13) relative some point away from the apex rather than the point 
Such point will ordinary point the differential equations, and 
there willbe loss generality itis chosen the value atrans- 
verse boundary the shell. order achieve this, let the distance 
and the length bah, are dimensionless parameters with 
and introduce new independent variable such that 


the differential equations (Eq. 13) now become 


2(m2 


+ 


«Theory Thin Elastic Shells,” 236. 
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and 


which the primes denote differentiation with respect and X(y), 
are obtained from Eq. making the substitution Eq. 15. The trans- 


Solutions Eq. relative the ordinary point are now sought 
the form 
Reasoning analagous that indicated the foregoing shows that solutions 
this form will converge for all <1, which includes the region interest 
for all conical shells finite length when the parameters and are proper- 
chosen. 
assumed that the surface loads can the form 


substituting Eq. and Eq. into Eq. 17, the following recurrence rela- 
tions are obtained: 
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for all 0,1,2,... Quantities with negative subscripts are taken zero. 

Eight linearly independent solutions the homogeneous system Eq. 

the coefficients u,, v,, from Eq. 20. The solution obtained from Eq. 

with these coefficients will denoted Similarly, 

can constructed setting all the leading coefficients equal 

zero except This process can continued until the solution Ug(y), 

has been obtained. particular solution Ug(y), can 

found the same way setting all the leading coefficients 

values. 

Hence, the general solution Eq. can written the form 


and 

The values ug, can adjusted satisfy system boundary 
conditions and When the values the leading coefficients ug, 
have been determined, then the solutions U(y), V(y), and W(y) satis- 
fying the boundary conditions can obtained from Eq. from Eq. 
making use Eq. 20. With this solution, the magnitudes the stress resul- 
tants are determined the following equations which are obtained making 
obvious substitutions into Eqs. and 
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NUMERICAL EXAMPLES 


Now there available least formal solution the indicated class 
problems, but the objective obtaining numerical results still some dis- 
tance away. The series Eq. converge slowly, and calculating the coef- 
ficients from the recurrence relations Eq. best slow and tedious 
task ordinary methods are used. order overcome these practical 
difficulties, IBM 650 computer has been used evaluate solutions. The 
calculation proceeds three stages which are outlined follows: 

given. Set and using the recurrence relations Eq. described 
previously, calculate the nine solutions etc., appearing 
Eq. 21. 

Stage II. Boundary boundary conditions must satisfied 
results Stage unnecessary use the computer develop state- 
ment conditions general, the boundary conditions lead eight 
non-homogeneous algebraic equations the eight unknown ug, V1, Wo, 
W1, which are solved using standard program for the computer. 

Stage III. Final Results.—Having determined the leading coefficients the 
power series solution (Eq. 18), from Stage II, these are now entered into the 
same program used Stage and the solution satisfying the boundary condi- 
tions evaluated for various values between The stress 
resultants are then calculated according Eq. 23. 

responds surface load 10-8 sin The shell assumed 
The results the calculation are given Table 1(a) and are 
plotted Fig. and Fig. 


corresponds surface load (3/ 16)(D/h3) sin The shell 
1(b) and are plotted Fig. and Fig. 

These examples are intended illustrate what can accomplished with 
the method described this paper. With the program which has been con- 
structed for the IBM 659, any meaningful values the parameters 
and can used, while the loads X(y), T(y), and Z(y) appearing can 
polynomials degree nineteen. should clear that problems with other 
solved the same manner. 

Calculations were carried out eight decimal digits and many 
terms were summed depending the value for whichthe series were be- 
ing evaluated. However, because the calculations involve subtracting quan- 
tities having similar magnitudes, our final results cannot correct eight 
significant figures. Independent calculations were made using different pro- 
grams for the computer and, least for the examples considered, the results 
were found consistent within 2%, 3%. far engineering appli- 
cations are concerned, apology needed for this degree accuracy since 
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FIG, 2.—SKETCH FOR EXAMPLE 
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FIG, 3.—RESULTS FOR EXAMPLE 
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the fundamental hypothesis which the theory based introduces inaccu- 
racy about this same magnitude. Because the expense operating the 
computer, the author was unable make more extensive numerical study 


the solutions, but this expected the object further research the 
subject. 
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ADDED MASS LENSES AND PARALLEL PLATES 


SYNOPSIS 


Contained herein study the added mass lens-shaped bodies, two 
parallel square plates, and two parallel, infinitely long rectangular plates var- 
ious thicknesses. The added masses were determinded from mass-frequency 
relationship obtained immersing the objects water and accelerating them 
oscillatory motion. The experimental results are good agreement with 
the analytical studies corresponding potential flows. 


INTRODUCTION 


The motion solid body through infinitely extended real fluid atrest 
is, general, accompanied the following three interdependent, well 
time dependent, occurences: boundary layer, wake, and fluid mass transport. 
Consequently, the object subjected resistance. steady motion, bound- 
ary layer and the characteristics the wake and, hence, the total resistance 
dependon the velocity, some physical properties the fluid, andon the geo- 
metrical characteristics the object. Although, the fluid mass transported 
the mass fluid entrained the body, including the wake-induced mass, also 
depends onthe cited variables, does not contribute totalresistance. un- 
steady motion, however, the total resistance affected several ways. The 
velocity dependent drag affected that not for steady flow. 
inertial force occurs because the acceleration the additional mass en- 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Assoc. Prof. the Engrg. Mechanics Dept., Univ. Nebraska, Lincoln, Nebr. 
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trained the body. Finally, the geometry the wake varies with time, and 
both the drag and the inertial force are affected markedly. For oscillatory 
well unidirectional, unsteady motion, therefore, the analytical determination 
the inertial force and the drag appears very difficult. Studies? made 
the vortex shedding seem throw some light the understanding the gen- 
eral problem. 

The inertial force can, however, computed for the irrotational flow pro- 
duced ideal incompressible fluid the vibratory motion immersed 
object. The usefulness the preceding theoretical prediction, inertial force 
the inertia coefficient, obviously limited the very severe restrictions 
made its evaluation. the moment zero velocity there wake and 
the flow pattern closely resembles that for irrotational flow without separation. 
Therefore, the inertia coefficient obtained flowtheory consti- 
tutes single pointon the general drag-inertia-coefficient curve,2 whichis the 
approximate starting point motion. Approximate, because the fluid accumu- 
latedin the boundary layer prior separation causes the actualinertial forces 
slightly higher than what the theory would predict. 

Since the time Buat, the added mass has been the subject many ex- 
perimental and analytical investigations. For three-dimensional flows how- 
ever, the added mass has been determined explicitly only for very limited 
category geometrically prescribed objects. For all elliptic cylinders with 
one axis parallel the direction flow, the virtual mass equal the mass 
displaced circular cylinder for which the diameter the transverse axis 
bouchinsky for cylinders with rectangular cross sections. Also available 
terms integrals are coefficients for any ellipsoid, provided only that 
one axis parallel the direction For all these cases, irro- 
tational flow without separation assumed. 

Because separation occurs flow past any body which not well stream- 
lined, the theoretical results are not necessarily useful analyzing flows past 
structures conventional design. more realistic flow past 
flat plate was Riabouchinsky, second plate downstream 
the first closes the flow. The cavity, bounded the two plates and the con- 
necting streamlines, mathematical approximation actual wakes. 
Birkhoff al. extended® this analysis and presented drag coefficients for var- 
ious cases, The method was again extended recently include the evaluation 
the virtual mass systems with various spacings the 
Other ways achieving finite cavities are the reentrant Gilbarg and 
Rock and the parallel after-body8 Roshko. 


“Vortex Formation and Resistance Periodic Motion,” McNown and 
Keulegan, January, 1959, pp. 1-6. 
“Hydrodynamics,” Lamb, Dover Publications, Inc., New York, 6th Edition, 
1945, 

Resistance des Fluides,” Riabouchinsky, Comptes Rendus, Conres 
International des Mathematiciens, Strasbourg, 1920, pp. 568-585. 

“Wall Effects Cavity Birkhoff, Plesset and Simmons, Quar- 
terly Applied Mathematics, January, 1952. 

“Note Estimation Virtual Mass Coefficient,” McNown, (April 17, 1955, 
Unpublished, Sandia Corp.) 

Naval Ordnance Lab. Memorandum 8717, Gilbarg and Rock, 1945. 


New Hodograph for Free-Streamline Theory,” Roshko, NACA TN-3168, 
1954, 
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The classical theory without separation appropriate for motions which 
the displacement small. has been used the forces ondams 
during earthquakes, 9,10 for studies bodies vibrating liquids, 11,12 and 
for computations gust and vibration loadings airplane wings. Considerable 
discrepancy found between the results this type and those obtained for 
flows which the amplitude the motion large enough for separation 
pocket occur. wake-induced mass which varies continuously during the 
period vortex shedding had taken into consideration. 

the present study, the experimentally determined added-mass coefficients 
for lenticular object formed two intersecting spheres, for two completely 
separated spheres, for two parallel square plates, and finally for two parallel, 
infinitely long rectangular plates various thicknesses have been added the 
existing list added mass coefficients. 

this time, appropriate discuss briefly the evolution the idea 
added mass. Definitions virtual mass have been presented from several 
points view all leading the same result. Lamb presented3 method 
which the integral the Euler equation, including the term for unsteady flow, 
used. Another method integration throughout the fluid the in- 
cremental force which the product elementary mass and the local ac- 
celeration. Then, the quotient the total force required produce the accel- 
erations throughout the fluid divided the acceleration the body defines the 
virtual mass. The usual derivation virtual mass, however, from kinetic 
energy which suffers from the inherent weakness that assumes constant 
velocity for the confined solid body, whereas all outward manifestations mass 
are associated with acceleration. Therefore, velocity must allowed vary 
and the fact should recognized that the added mass can yield momentum 
well energy. Sir Charles Darwin has that object moving 
infinitely extended ideal fluid, besides pushing the particles aside temporarily 
passing, also displaces the fluid particles permanently the direction its 
motion and that this displaced mass fluid enclosed between the initial and 
final positions fluid particles, fact, the added mass itself. Although, 
Darwin’s physical interpretation ofthe added mass quite what 
was known, does not any way disprove the well-known kinetic-energy con- 


cept nor does constitute simpler method determine the added-mass co- 
efficients theoretically. 


THEORY 


Since the time rate change kinetic energy equal the rate doing 
work, upon integrating, there results 


Pressures Dams During Earthquakes,” Westergard, Trans- 
actions, ASCE, 1933. 

“Hydrodynamic Earthquake Forces Submerged Structures,” McNown, 
Proceedings, Third Midwestern Conference Fluid Mechanics, 1953. 

“Virtual Mass and Acceleration Fluids,” Stelson and Mavis, 
Transactions, ASCE, Vol. 22, 1957, pp. 518-530. 

“Virtual Masses Rectangular Plates and Parallelpipeds Water, Yee-Tak 

Yu, Journal Applied Physics, November, 1945. 


“Notes Hydrodynamics,” Sir Charles Darwin, Proceedings, Camb. Phil. Soc. 
49, 1953, pp. 342-354. 
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For three-dimensional flows the velocity potential has been determined ex- 
plicitly only for very limited category geometrically prescribed objects. 
The added-mass coefficients two separated spheres andof the lenses are the 
only ones determined analytically among those herein. The flow 
studied further. The added-mass coefficient two parallel, separated, 
square, and infinitely long rectangular plates have not yet been determined an- 
alytically. The latter problem should not confused with that one treated 

classical procedure for finding the flow due the vertical motion the 
two spheres bythe method images. The flow produced single sphere 
moving alone infinite body fluid that created dipole suitable 
strength placed the center the sphere. two dipoles are situated the 
centers two spheres, correct for their mutualinterference andto 
correct the destroyed boundary condition, one has correspond- 
ing images each dipole and continue this process untila convergent expres- 
the potential obtained. this potential functionis introduced 
Eq. and necessary integrations are performed, one for the added 
mass coefficient. Throughout this study Cis referred the displaced mass 
the body. 

For two completely separated spheres, the added-mass coefficient obtained 


inthe manner described foregoing deduced fromthe Shiffman 
and Spencer: 


When the spheres are tangent, 


The added mass for lens-shaped body between two intersecting identical 
spheres given by17 


which 


“Journal fur Mathematics,” Mehler, 68, 1868, pp. 134-150. 
ZaMM 12, Wagner, 1932, pp. 193-215. 
E., Tech. Note Aero., Monaghan, 1989, 1944. 
“Flow and Incompressible Fluid about Lens,” Shiffman and Spencer, 
Quarterly Applied Mathematics, Vol. No. 1947, pp. 270-288. 
arterly Applied Mathematics, Payne, Vol. 1952, pp. 197-204. 
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for odd, and 


for even and for 


The results obtained from these and from previous equations have been re- 
calculated and presented Fig. together with the rate change the added- 
mass coefficient with respect tothe ratio B/R. Thetheoretical results obtained 
from Riabouchinsky’s for two dimensional rectangular prisms are in- 
dicated dashed lines Fig. 


EXPERIMENTAL PROCEDURE 


The experimental apparatus was designed impart vibratory motionto con- 
nected pairs spheres, parallel square and rectangular plates, and lens- 
shaped bodies such manner that all stages motion the kinematic and 
kinetic characteristics the experimental system and test objects could easily 
deter mined. 

picture the apparatus one sample object from each group tested 
are shown Fig. The test bodies were immersed rectangular tank, 
in, square section, and in. deep, filled with tap water. in. long small 
aluminum beam was simply supported pivot bearings very carefullyto elim- 
inate end restraints. the mid point and below the beam elliptical alumi- 
num ring (major axis in., minor axis in., thickness 0.10 in., width 0.50 in.) 
was securely attached the beam with the major axis the horizontal posi- 
tion. This ring constituted the elastic element for measuring directly. 
indicate the ring deformations two pairs SR-4 strain gages were glued 
the ends the major axis, the inside and outside the ring. The four 
strain-gage the bridge which led universalanalyzer. 
having four strain gages the ring the sensitivity was increased and cor- 
rections were needed for temperature changes. analyzer is, 
turn, relayed one the channels two-channel magnetic oscillograph. 
The ring was calibrated statically for both tension and compression before 
mounting the beam. Thereafter, the calibration was repeated before and 
after each run the mountedposition make sure that accidental changes 
occurred the strain gages. The calibrations always reproduced and plotted 
straight line within the loading range. Furthermore, the system was cal- 
ibrated dynamically. known mass suspending air was attached rigidly 
the ring and then the beam was vibrated means magnet-coil system. 
The frequencies were measured with sonometer and, also, with universal 
counter Accelerations, displacements, and velocities were meas- 
ured with vibration meter (model Separate measurements the latter 
kinematic variables also providedanother check inthe measurements since the 
velocity, displacement, acceleration, and frequency are related with simple 
well-known expression harmonic motion. Using the force recorded the 
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oscillograph and the acceleration measured, the total mass the suspending 
rodand the mass the sample body was calculated. the sus- 
pending rod was previously determined and remained constant throughout the 
experiment, the mass the sample object was determined and compared with 
the actual mass determined through the measurement its weight. Such 
comparison always gave values within the actual mass and was regarded 
satisfactory view the complexity the system. The frequencies se- 
lected were close neither fundamental nor harmonic frequencies the 
beam nor the natural frequency the elliptical ring. The fundamental fre- 
quency the ring was computed approximately 850 cps. During each ex- 
periment oscilloscope tracings were carefully watched and was made sure 


that there were outside interference that would cause undesirable parasitic 
vibrations. 


3.—VIEW EXPERIMENTAL APPARATUS AND 
SOME THE TEST BODIES 


During dynamic calibrations was found that the air resistance entirely 
negligible. Following the calibration, the tank was filled with tap water and the 
object was vibrated the same frequencies and the accelerations were meas- 
From the force and acceleration readings the virtual mass was deter- 
mined. Damping caused water, which was quite small, included the 
computation the virtual mass. the tests, the velocities varied from 0.20 
in. per sec 0.60 in. per sec, the frequencies from cps cps, and the 
accelerations from in. per 400 in. per Finally, should 
noted that the experiments were repeated several each test object for 
various frequencies within the frequency range tested. For given object, 
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appreciable difference has been noted between the added masses determined 
using different frequencies, that is, from cps cps. 

Test objects were constructed brass, aluminum, and plaster paris. 
Square prisms consisted carefully sealed hollow aluminum blocks. Spheres 
were made brass shell, the inside which was filled with parafin wax. 
Lenses were made plaster paris and coated with parafin wax. The sizes 
the objects tested are noted Table 


TABLE 1.—DIMENSIONS TEST OBJECTS 


Object Dimension, inches 


Lenses 1.20 1.50 3.00 2.40 1.50 

3.00 3.00 4.00 2.00 1.00 
Square 


Parallel 


square 0.20 0.40 0.60 0.80 1.00 2.00 

prisms, 2.00 2.00 2.00 2.00 2.00 2.00 
szt For each, s/a varied with the connecting rod 

Parallel 

two- 

dimensional and same for parallel square prisms. 

rectangular Length in. 


prisms 


The added-mass coefficients for lenses and separated spheres obtained ex- 
perimentally and theoretically are shown Fig. Also illustrated Fig. 
arethe added-mass coefficients for square were obtained 
imentally only. The points representing the latter coefficients comprise the 
dot-dash line, shown Fig. This dot-dash line used later Fig. 
define the asymptotic curves for various values the parameter t/a. Also 
illustrated Fig. plot the rate change the added-mass coeffi- 
cient for lenses with respect B/R. first, noted that the experimental 
values are satisfactory agreement with the theoretical values. The second 
point notedis that the added-mass coefficient two intersecting spheres 
minimum about B/R 1.75. one can consider two such intersecting 
spheres the equivalent potential model sphere and its wake together 
with two symmetrical vortices situated the wake, one finds that during the 
unsteady motion the added-mass coefficient the model drops from 0.50 
0.30. not the writer’s intention replace the complicated unsteady mo- 
merely point out the fact that the added-mass coefficient drops considerably 
the wake grows longer. However, immediately shedding the vor- 


tex, the motion becomes more complicated and present nothing much can 
said about the variation 
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Figs. and curves with dotted lines ratio added mass dis- 
placed mass for unslotted, two dimensional, right rectangular prisms and for 
full right square prisms, respectively, function s/a, for constant values 
t/a. Curves with full lines give the aforementioned ratio for the separated 
and rectangular plates function the same dimensionless ratios. 
the dotted lines Fig. are obtained from Riabouchinsky’s 
and the dot-dash lines Fig. are derived from the experimental 
curve shown Fig. The direction the vibratory motion also shown 
the figures. 

examination the figures reveals that the plates are separated from 
each other but very little distance, the added-mass coefficient decreases 
slightly, that is, the fluid between the plates behaves were part the 
solid full parallelpiped. the plates are set farther apart, the added-mass 
coefficient increases, but always remains smaller than the added-mass coef- 
single plate. Hence, when the distance between the plates not too 
large comparison with the width the plate, there mutual interaction 
the flow patterns around the two plates. the plates become farther apart, 
action becomes negligible and each plate behaves independently. 

should noted that although the basic concepts the subject matter 
contained herein have been expressedin similar ways many others, the ex- 
perimental for the added mass multiple spheres, lens-shaped 
bodies, and parallel square and rectangular plates new contribution. 
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NOMENCLATURE 


half the thickness lens measured the direction motion, 
(Fig. 1); 


added mass coefficient, (added mass/displaced mass); 
function; 


displaced mass; 


radius the circular plane surface generated the intersection 
two spheres. (Fig. 1); 


radius sphere; 
see Fig. 


kinetic energy; 
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velocity the undisturbed stream; 


coordinate axes; 


half the exterior angle between the two intersecting spheres 


density fluid; and 


velocity potential function 
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STRESSES DUE THERMAL GRADIENTS REACTOR SHIELDINGS 


SYNOPSIS 


The stresses due sudden rise surface temperature one side 
fixed ended arch with rectangular cross-section are evaluated, taking into ac- 
count the temperature variation the arch. The results are ex- 
tended the case surface temperature which first increases linearly and 
then remains constant, means Duhamel integral the case arbi- 
trary temperature applied the surface the structure. 


INTRODUCTION 


Nuclear reactors require shielding structures capable serving their 
structural integrity environment which differs many respects from 
that other civil engineering structures. The necessarily massive dimen- 
sions the shielding make particularly sensitive temperature changes 
andits structural design often primarily governed thermal stresses rather 
than dead load and live loads. Thermal stresses can arise from: 


Slow, cyclic, temperature variations due natural causes. 
Slow temperature variations due normal operations. 
Slow heat generation due absorption gamma rays. 


Note.—Discussion open until November 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Chf. Engr., Paul Weidlinger, Cons. Engr., New York, Y., and Adjunct Prof. 
Engrg., Columbia Univ., New York, 

Assoc., Paul Weidlinger, Cons. Engr., New York, Y., and Civ. Engrg., 
Columbia Univ., New York, 
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rises (thermal shock) due accidental causes the 
reactor operation. 


Temperature variations associated with and are normally sufficiently 
slow allow steady state conditions with corresponding linear gradient within 
the shielding assumed. The determination the state stresses can 
thus obtained well known engineering methods. 

Case leads nonlinear steady-state temperature gradient 
determined the energy absorption coefficient the material known. The 
corresponding state stress canbe found the equations given the section 
strains due thermal stresses. 

somewhat more difficult case because the nonlinear 
transient temperature gradient due the thermal shock cannot given 
closed form valid for any time interval. the other hand, the stresses the 
shielding structure due the kind high intensity shock considered here, 
rise very rapidly, and the interest the designer restricted the determi- 
nation the time when the stresses reach their greatest allowable value. Be- 
fore this time reached, provision must made reduce the temperature 
within the shielding cooling, since most instances would uneconomi- 
cal design the shielding sustain very high intensity temperatures for ex- 
tended periods. Since the maximum stress early time tensile stress 
that strengthof the concrete not exceeded beyonda depth 
equal the coverage over the reinforcing bars. Under these conditions, the 
surface the shielding may crack spall, but its essential structural in- 
tegrity preserved. 

For convenience, the time-temperature history the inside surface the 
shielding first considered step function, corresponding thermal 
pulse zero rise time, and the corresponding temperature distribution the 
arch evaluated. However, the containment sphere normally separated 
from shielding air other insulation, the temperature 
pulse will have finite rise time with linear non-linear variation during 
this rise. Beyond the rise time, the surface temperature reaches and main- 
tains maximum value the time period interest. 

Once the temperature distribution the arch has been evaluated for the 
case appliedstep pulse (zero rise time), these results can usedas in- 
fluence coefficients determine the corresponding temperature distribution 
arch for appliedarbitrary temperature its inside surface, means 
Duhamel integrals. Three cases applied temperatures are considered 
the subsequent sections: step pulse with zero rise time; pulse with 
linear rise time; and pulse with arbitrary nonlinear rise time. 

The results obtained the present study apply the supporting arches 
reinforced concrete ribbed dome subjected the above described tempera- 
ture impacts their inner surface. The outer surface remains constant 
temperature during the time period interest. 

Arches shallow depthare considered, where the influence the curvature 
the temperature gradient can Direct due 
the nonlinear gradient, together with end rotations and displacements the 
arch are determined. Using these results, the end moments and thrust the 
arch can obtained for any support condition full partial fixity. nu- 
merical example for built-in symmetrical arch presented. 


Stresses Curved Beams,” Bruno Boley and Barrekette, Journal 
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REACTOR SHIELDINGS 
TEMPERATURE GRADIENT DUE THERMAL SHOCK 


Step Function with Zero Rise Time Applied Inner Face free 
ended rectangular beam depth 2c, and length (Fig. 1), insulated its 
lateral and end boundaries, has its boundary suddenly raised tem- 


following 


Short Time Solution.— 


Complete 


(y, 


sin 


which (in square feet per hour) the thermal diffusivity the material 
the beam and the function erfc the error function complement 


Eq. shows that (y, asymptotically linear distribution. For 
early times the series Eq. converges rapidly that for most practical 
cases, well represented the first term, 

The temperature distribution Eqs. can used with good approxi- 
mation for arches, provided their radius curvature large compared 
their depths 2c.3 

typical case, Fig. shows the temperature distribution (y, t)for 
reinforced concrete archwith rectangular cross section, 5-ft deep. The value 
the thermal diffusity “k” was taken 0.33 per hr. 

After 16-hr the maximum difference between the steady state temperature 
and (y, t), which occurs less than 20%; after hr, less 
than 2.5%. The error when only the first term Eqs. and 
used, less than for hr. 

The temperature distribution actual beam arch, which not insu- 
lated the ends, not function the coordinate only, since there heat 


Not the neighborhood the ends the beam. 


“Conduction Heat Solids,” Carslaw and J.C. Jaeger, Clarenden Press, 
1947, 250 ff. 
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flow the ends the axial direction Hence the temperature distribution 
except the neighborhood the ends. This neighborhood small 2c/R 
small, which the case for “long” arches. 

Arbitrary Temperature Pulse T(t) Applied the Inner Face the Arch.— 
Let t)be the temperature distribution the arch due unit step tem- 
perature pulse applied the inner face the arch, that is, let given 
Eq.1 with set equal unity. Using these values influence coefficients, 
the temperature distribution T(y, produced arbitrary temperature pulse 


FIG, 1.—GEOMETRY 


T(t) applied the inner face the arch can evaluated means Du- 
hamel integral: 


evaluated analytically, numerically. general, only the first term 

the series expansion for need used, that is, 


Depending the complexity the function the integral Eq. can 


thus considerably simplifying the integral Eq. 4b. 
Temperature Pulse T(t) With Linear Rise Time Applied the Inner Face 
the Arch.—For the particular case linear pulse, 
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where constant temperature and particular value (Eq. 11), 

Eq. may integrated analytically term term using the relation 

vt-T 


For practical purposes, only the term need kept the expression 
for (y, t), Eq. becomes 


2 vk (t - 7) (t - 7) 


upon integration, 


7 oO 1 1 


+ (y) * to fo erfc fo (y) 


which 


advantageous some problems practical interest consider the 
distribution (y, due the pulse shown Fig. 


the function (v. obtained from the relation 
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STRAINS DUE THERMAL GRADIENTS (y, 


has been shown Boley and that the longitudinal strains 
curved rectangular beam radius and depth can well approxi- 
mated the strains straight beam, for cases which the depth radius 
ratio the curved beam small. The error the strain that arises from 
neglecting the effect the curvature less than for 0.2. For practi- 
cal applications, many reactor shielding structures fall 
the effect the curvature the beam will neglected the following de- 
velopment. 


T(t) 


t/te 


FIG, 3.—TEMPERATURE 


The longitudinal strains due temperature distribution (y, t), ina 
beam rectangular cross section given the 


where (with denoting the width the rectangular cross section) 


(cross section area beam) (14) 


and (in feet per degree the coefficient thermal expansion 
the material the curved beam. The last term Eq. the strain which 
produces stress, due the thermal expansion, thefiber adistance from 
the neutral axis. The first three terms are the strains due restraint 
cised the adjacent fibers the y-fiber and which produce the elastic 
EX. The integrals Eq. are extended over the entire cross sectional 
area the beam, leading the expression, 


Determination Temperature Stresses and Deflections Two-Dimensional 
Thermoelastic Problems,” Bruno Boley, Journal the Aeronautical Sciences, Vol. 
23, No. January, 1956. 

Elasticity,” Timoshenko, McGraw-Hill Book Co., New York, 1954, 
204. 
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The the upper and lower fibers, respectively the arch may 
evaluated and will used the next section todetermine the endrotations 
and displacements the arch when the ends are considered free. 


END ROTATIONS AND DISPLACEMENTS 


Noting that the temperature the upper face the arch kept zero 
degrees, the strains the upper and lower fibers, and respective- 
ly, become, using Eq. 15, 


(16) 
where 


and the maximum value which the temperature pulse main- 
tained. 


(Fig. 4), are given 


Since the formulas for strain, Eqs. and 15, have been derived the as- 
sumption that sections plane before heating remain plane after the 
rotation given section along the arch given 


the end rotation becomes: 


Under the assumption small displacements, the slope the thermal de- 
flection curve z(x) (Fig. equals that 
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Ec 


and integration, the end displacement right angles the arch, given 


The corresponding horizontal displacement the end the arch central 
angle becomes (Fig. 


The rotation and displacement “a” have been derived considering the 
arch unrestrained its ends. End moments and horizontal shears must 


FIG, 4.—GEOMETRY ARCH 


now applied the ends the arch force them back into their original 
fixed position. 


END MOMENTS AND SHEARS 


Let and the moment and shear required produce unit rotation 
and zero horizontal displacement the end arch length 2L. Also, let 
and the moment and shear required produce unit horizontal dis- 
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placement and rotation the ends the arch. The four quantities hg, 
and can easily determined for any given arch any the usual 
structural methods available the designer. The moments, and shears, 
required produce end rotations and displacements, equal and opposite 
those produced the thermal gradient and a), are given the equations 


M+h, H+a=0 
(hg 


Eqs. are statement that the ends the arch are fixed. similar meth- 


ods, partially fixedor hinged arch moments and shears can obtained. From 
Eq. 25, 


2 hg a- hg 6 (26) 
Mg hy - hg 


FINAL STATE STRESS THE ARCH DUE THERMAL SHOCK 


The final state stress the arch, produced the temperature pulse 
tothe inner face -c, obtained the superposition (a), the 
state stress corresponding the strains Eq. 13, and (b), the state 
stress produced the moment (Eq. 26) and the shear (Eq. 27). The state 


where the modulus elasticity the arch. may written 
terms the quantities and 


RESULTS FOR STEP PULSE, T(t) 


Substituting the term Eq. into Eqs. and 18, expressions for 
and are derived terms the integrals the error function complements® 
(ierfc): 


Stresses Curved Beams,” Appendix “The Error Function and Re- 
lated Functions,” Bruno Boley and Barrekette, Journal Aeronautical Sciences, 
Vol. 25, No. 10, Oct., 1958. 
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FIG, 5.—STRESS DISTRIBUTION 
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and 


Once and have been obtained for given geometry, the values 
and “a,” given Eqs. and respectively can computed; hence, and 
from and 27. 

The total stress the arch, (y, t), due the temperature step pulse 
the inner face, can obtained superposition (a), the direct thermal 
stresses, calculated from Eqs. and and (b), the direct stress due the 
application the moment and the shear the ends the arch. 

For purposes illustration, the direct thermal stresses, computed from 
Eqs. and are shown for the archof rectangular cross section (depth ft, 
0.33 per hr.) which was discussed under the section arbitrary 
temperature pulse applied the inner face the arch. Fig. gives the direct 
1/6hr, 1/4 hr, and 1/2 hr. 

Analogous results can obtained for other temperature variations com- 
puting the values and (Eqs. and numerical integrationonce the 
value T(y, has been computed from 4(a), 4(b) For temperature 
inputs practical interests, the integrands are reasonably well behaved and 


standard numerical integration techniques such Simpson’s Rule? may 
employed. 
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STRESSES AROUND RECTANGULAR OPENINGS PLATE 


Boyd Phillips,! and Ira Allen? 


SYNOPSIS 


Curves have been developed for determination the normal stress around 
rectangular openings when subected biaxial stress fields. Height 
width ratios and 2.5 have been studied. The study was made ex- 
perimentally utilizing the photoelastic interferometer method stress analysis. 


INTRODUCTION 


the design massive structures, suchas mass concrete dams, the prob- 
lemof determining the state stress around rectangular openings frequently 
encountered. The opening may vertical elevator access shaft, hori- 
zontal sloping gallery, other such opening the mass. 

The stress distribution around openings uniform stress field has been 
investigated However, those studies werefor uniaxial stress 


Note.—Discussion open until November 1960. extend the closingdate one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechaincs Division, Proceedings the 
American Society Civil Engineers, Vol. 86, No. June, 1960. 

Engr., Design Div., Bur. Reclamation, Denver, Colo. 

Engr., Design Div., Bur. Reclamation, Denver, Colo. 

Totale Zugkraft and Offnungen Einem Einachsigen Druckspannungsfeld,” 


Hiltscher and Pant, Der Bauingenieur, Vol. 32, No. 12, December, 1957, pp. 
470-474, 


«Stresses Around Gallery-Determined the Photoelastic Interferometer,” 


Phillips and Zangar, Proceedings, Soc. for Experimental Stress Analysis, Vol. 
VIII, No. 1951, pp. 187-208. 


“Stress Analysis Applied Underground Mining Problems—Part I,” 
Dept. the Interior, Bur. Mines, Report 4192, March, 1948. 
«On the Stresses the Neighbourhood aCircular Hole Strip Under Tension,” 


Howland, Philosophical Transactions, Royal Soc. London, Sec. Vol. 
229, 1930. 165 
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field. This study was undertaken simplify the determination stress dis- 
tribution for the design reinforcement for rectangular openings various 
height width ratios and for combination stress fields. 

For the preliminary design reinforcement around gallery other 
opening mass concrete dam, usually necessary know only the net 
tensile force acting onthe center lines the opening. the caseof horizon- 
tal gallery the tensile force caused the predominant dead loadof the struc- 
ture, and may reduced horizontal compressive stress field caused the 
reservoir load the arch load the case arch dam. For the final de- 
sign, stress determination along various lines radiating fromthe opening due 
the two principal stress fields may required. 

Notation.—The letter symbols adopted for are defined where 


they firstappear andare arranged alphabetically, for convenience reference, 
the Appendix. 


APPLICATION DATA 


Fig. and are curves for normal stress perpendicular the center 
line rectangular openings (having ratios 0.4 2.5) infinite 
plate due uniform stress fields parallel perpendicular the line. 
Fig. and indicates the intensity the stress field acting the opening 
the direction indicated the subscript. Figs. and are photoelastic in- 
terferometer stress analysis with B/A from 0.4 2.5 and D/A from 0.0 
0.42. Fig. shows the stress field with along the Y-axis and Fig. shows 
the field with along the Y-axis. was found inthis studythat the tensile 
stress normal tothe center line dueto stress field parallel tothe center line 
increased slightly with increase the B/A ratio. The increase was 
slight, however, that for design purposes the same stress curve may used 
for all B/Aratios studied (Fig. 1). Other investigators3 have also 
same conclusion. The compressive stress the center line the boundary 
the opening due stress field perpendicular decreases rapidly with 
decrease B/A, falling below the maximum B/A ratio approxi- 
mately 0.5 (Fig. 2). 

find the net tensile force acting the center line rectangular open- 
ing due biaxial stress field, use the following procedure: 


Determine B/A and 

Using this B/A value, tabulate the values from the curves Fig. 
Multiply the values 


Determine the tensile area between the two curves. 


This the net tensile force terms and 


Asa numerical example, consider horizontal rectangular gallery mas- 
sive structure having gallery width ft, height 1/2 ft, vertical com- 
pressive stress field 100 psi, and horizontal compressive stress field 
psi. The solution this example shown Fig. 

Since rectangular openings having B/A ratio 1.5 are common mass 
concrete structures, extensive analysis was made the stresses around 
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PLATE OPENINGS 


For design purposes the normal stress on the 
centerline parallel to the stress field direction 
may be considered to be independent of B 

For p, = 0, the tensile area = 0.149 Ap, 

0, is the normal stress perpendicular to Y-axis 

D is the distance from the boundary of the 
opening along the Y-axis 


TENSILE STRESS 


0.24 
D/A 


FIG, 1.—RECTANGULAR OPENINGS INFINITE PLATE 
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a) opening along the Y-axis. 
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openings this ratio. The results this phase are shownon Fig. 
Table 

other methods analysis the normal and shear stresses, 
are determined’ the point the structure where the opening tobe located. 
They areassumed stress fields their respective intensities and 


Notes: 0, is the normal stress perpendicular to Y-axis. 
4 D is the distance from the boundary of the 1.9 


B/A 


FIG, 2.—RECTANGULAR OPENINGS INFINITE PLATE 


directions, px, Py, and Sxy, acting the opening. determine the normal 
stress, the lines around the opening, multiply the stress coefficients tab- 
ulated Fig. their respective intensities and sum the three resulting 


Load Method Analyzing Arch Dams,” Boulder Canyon Project Final Re- 
ports, Part Tech. Investigations, Bulletin Bur. Reclamation, 1938. 
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A = 5 feet 

B = 7.5 feet 
B/A 1.5 

py = 100 psi 

p, = 10 pst 

k = /P, = 


TABULATION FROM CURVES 


ON FIGURE 2 

0.02 0.177 

0.166 
0.8 


The shaded area between 
the two curves 
0.0912 Apy 
0.0912 (5) (100) (144) 
6570 Lb/ft-length gallery 


The net tensile force 


--- Tensile stress, On /Py due to vertical stress field 
(From Figure 


ie) 0.04 0.08 0.12 0.16 0.20 0.24 0.28 O.32 0.36 0.40 
D/A 


‘Compressive stress, kO;, / p, , due to horizontal 
stress field. (From above tabulation) 


FIG, 3.—USE CURVES FIGS. AND 
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COEFFICIENTS 


Notes: D is the distance from the boundary of the opening along the line. 
0, is the normal stress perpendicular to the line. 
compression. 
is tension. 
On/Py and O,/p, stresses are symmetrical about both centerlines. 
stress signs are shown for and 
ond opposite for IT and IV. e 


FIG, 4.—PHOTOELASTIC INTERFEROMETER STRESS ANALYSIS FOR B/A 1.5 
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stress effects. instead ox, Oy, and Txy (the shear stress acting Fig. 
the principal stresses, and and their direction, are given, and 
may rotated into and Txy. 

anumerical example the use the stress coefficients, consider rec- 
structure with maximum compressive principal stress field 100 psi 
1/2 counterclockwise from vertical, anda minimum compressive principal 
stress field psi. 

Given then that; ft, 1/2 ft, B/A 1.5, psi, 100 psi 

Using the rotation equations: 


obtain the following values: 23.2 psi, 86.8 psi and, Txy 
For the tabulation see Table 


TECHNICAL DETAILS 


The complete state stress any point two-dimensional problem 
defined when the two principal stresses and their directions are known. The 
photoelastic polariscope provides only and requiring another meth- 
od, such graphical integration, lateral extensometer, membrane analogy, 
electrical analogy tray, determine free boundary where one 
principal stress zero and the direction the other parallelto the bound- 
ary, the polariscope provides very rapid method determining boundary 
stress. For determining interior stresses, however, the photoelastic inter- 
ferometer® much more convenient because provides and directly. 
addition, less load required obtained with the po- 
lariscope, and time-edge effect has appreciable influence the optical meas- 
urements. 

The for the model was allyl diglycol carbonate resin. The 
variable height the opening was from 0.2 in. 1.25 in. with constant width 
was centered 5.8- 5.8-in. plate 0.25-in.-thick. The open- 
ing with the B/A ratio 1.5 was 0.6 in. 0.9 in. the same size plate. 
compressive load approximately 500 psi was applied the model through 
pair multiple-element fluid load shoes, developed recently the authors. 
Fig. photograph the model and the loading mechanism. 

effort evaluate the relation the normal stress around opening 
plate finite dimensions that infinite plate, the normal stress 

long the center lines parallel and perpendicular the stress field direction 


«photoelastic and Experimental Analog Procedures,” Moody and 
Phillips, the Interior, Bur. Reclamation, Engrg. Monograph No. 23, 
August, 1955. 
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were determined for series circular openings 5.8-in. plate. 
analytical solution exists for the case circular opening infinite 
center line parallel tothe stress field normal stress 
experimentally was the analytical value until the diameter 


1.—RECTANGULAR OPENINGS INFINITE PLATE 


Reverse sign for 
In psi. 


the opening approached in., beyond which the experimental value exceeds 
the analytical value. the center line perpendicular tothe stress field direc- 


Elasticity,” Timoshenko, McGraw-Hill Book Co., New York, 
1934, 


a 
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tion, the normal stress increased near the opening with increasing opening di- 
ameter manner that could approximately related the average normal 
the opening thanthe analytical curve, falling below edge. 
The analytical curve, course, approaches asymptotically. some 
point then the two curves cross, and this point the correction for the finite 


FIG. 5.—RECTANGULAR OPENINGS INFINITE PLATE 


plate effect would zero. This point varies with opening diameter, but for the 
corrections applied the lines perpendicular the stress field direction 
Half the correction was applied the lines the direction the stress 
field and correction the lines parallel the stress field direction. 
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APPENDIX. —NOTATION 


Coordinate System Positive Sign Convention 
FIG. 6.—NOTATION SYSTEM 


Normal component stress parallel the X-axis acting with 
opening present. 

Normal component stress parallel the Y-axis acting with 
opening present. 

Shear stress acting with opening present. 


Normal stress perpendicular the line being considered. 


Principal stresses acting with opening present. 

Angle the principal stress measured counterclockwise from the 

Intensity the stress field acting the opening the X-direction. 
Intensity the stress field acting the opening the Y-direction. 
Intensity the shear stress field acting the opening. 

Intensity the stress field acting the opening the direction 


Intensity the stress field acting the opening the direction 


Dimension the opening the X-direction. 
Dimension the opening the Y-direction. 


Distance from the boundary the opening along the line being con- 
sidered. 
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ELECTRICAL ANALOG FOR EARTHQUAKE YIELD 


Discussion John Blume 


JOHN ASCE.—The research work reported Messrs. By- 
croft, Murphy, and Brown particular interest the writer because his 
having worked the same problem from different approach. The efforts 
these investigators with the electrical analog constitute definite contributions 
the storehouse available informationin earthquake research, particularly 
inelastic dynamic response about which relatively little has been published. 

There can little remaining doubt that the history damage, and lack 
damage, severe earthquakes cannot reconciled unless energy loss and 
inelastic behavior are included major factors.2 Viscous damping, rubbing 
friction, and equivalent viscous composite materials and building 
elements are apparently insufficient magnitude and effect reduce elastic 
spectral responses from actual severe earthquakes the yield values even 
the ultimate strength values buildings.4 This true for modern buildings 
designed under seismic codes, well for the traditional type buildings. 

much more logical procedure than design elasticity alone 
design, refinement design, inelastic procedures that take into con- 
sideration the energy absorbing value all materials and building elements. 
simple method for doing this normal structural design terminology has 
been presented the writer.4 energy technique,” general, 
thatit providesfor any andall combinations ductile, brittle, elasto-plastic, 
failing, work hardening materials and building elements, will applied 
herein the elasto-plastic system the authors used their research. 

The first step the procedure draw static load-deflection curve 
failure for the story building under consideration. should noted that 
such cannot plotted, any mathematical design procedure, matter how 
rigorous, merely rough estimate that may lead uneconomical designs 
inconsistent and unsatisfactory results actual earthquakes. The ad- 
ditional value most types stress and most building materials under 
rapid dynamic loading compared static loading roughly compensated 
for physically the gradual but declining rate deterioration resistance 
under the repeated cycles severe earthquake “loading” and yield excursions. 


Pres., John Blume and Assocs., Research Div., and Francisco, Calif. 

Dynamics Earthquake-Resistant Design,” John Blume, Pro- 
ceedings, ASCE, Vol. 84, No. July, 1958. 

“Frictional Effects Composite Structures Subjected Earthquake Vibrations,” 
Jacobsen, Dept. Mech. Engrg., Stanford Univ., March 1959, (California 
State Div. Architecture, sponsor.) 

Blume, Proceedings, ASCE, Vol. 85, No. September, 1959. 
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Thus, static load tests can reasonably used basis for constructing the 
required load-deflection diagrams that are actual approximations the spring 
factor the dynamic system. The spring factor may linear, hardening, 
softening, elasto-plastic, any combination. should plotted the point 
failure for the assumed initial force application. 

The “reserve energy reduction coefficient,” the basis reconciling 
any elastic design coefficient, with the elastic spectral earthquake response 
structure elastic period Stated another way, the factor that re- 
duces elastic spectral response accelerationto equivalent value for elastic 
design modified (by for more realistic inelastic response values. based 
upon equating the kinetic energy maximum response the total available 


strain energy and energy dissipation value. The for for single 
mass vibrating system is: 


the subscript referring the number amount any trial de- 
flection the structure unit under consideration, indicates the funda- 
mental period, elastic state, refers the adjusted period Aj, inelastic 
state, the elastic strain energy tothe design unit stress, the areaunder 
the diagram Ap. (Fig. 1), denotes the total area under the diagram 
trial deflection, and refers the “hump” area above straight line from 
the origin any point the diagram any deflection point, 

For the special (elasto-plastic) case hand, will considered equal 
Tg, since the elasto-plastic system does not materially change its spring 
factor. The “hump” deduction, design allowance for the possible soften- 
ing deterioration resistance under repeated cycles), will considered 
zero for this assumed perfectly elasto-plastic with deterioration. 
Thus, the expression becomes simply 


These values are initially assumed general for any earthquake spec- 
trum and for any value damping. This will evaluated herein. The only 
remaining step toconsider the specific earthquake spectrum for any select- 
damping factor. Letting equal the elastic ratio 
the design coefficient, the structure barely adequate energy dissipa- 
tion value (at safety factor), deflection not necessary assume 
value only the yield point, since with more permanent yield set (greater 
ductility factors) the structure can absorb much more energy. Thus, the au- 
thors’ Fig. does not indicate all important parameters, since the equivalent 
elasto-plastic accelerations would vary with the deflection beyond the yield 
the energy dissipation capacity. was decided plot ductility versus 

Fig. was utilized obtain data for computing the authors’ ductility fac- 
tors, deflection yield-point units. The maximum yield displacement was 
added the elastic displacement, and the sum divided thus 
obtained for 0.20 are shown the first four columns Table The fifth 
and sixth columns show the determination the equivalent “R” value from 
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Fig. 16, and the assumption that the design unit stress, including one-third 
increase over normal (for earthquake loading), was 0.75 times the yield unit 
stress, which reasonable assumption for “elasto-plastic” mild steel. 
Columns and Table showthe same data asfor columns and but for 
another source elastic data furnished the writer the authors part 
another research The variation the values columns and 
noted. Columns and are plotted Fig. together with other data. 

Fig. plot elasto-plastic ductility factor versus the ratio elasto- 
plastic spectral elastic spectral acceleration, the “reserve 
energy reduction coefficient.” The heavy solid line Rj, calculated, pre- 
viously described, the reserve energy procedure and Eq. noted that 
the decrease from unity design value the yield point due the 
safety factor lineal basis, and beyondthe yieldpoint due 
energy dissipation. Fig. indicates these variations the elastic response 


TABLE 1.—DETERMINATION AND DUCTILITY FACTOR 


Period, 
sec. in. Yield Acceler- 0.75 0.24 
ment, factor Ratio 
in, 


(3) (5) 


and right-hand columns Table show for the two sets 
analog data, both which are plotted Fig. for the corresponding ductility 
factors. The general relationship tothe curve fairly close for the authors’ 
values, although the vertical scale variation between the new research pro- 
previously referred and Fig. data considerable. 

The discussion the authors’ paper Messrs. Berg and Thomaides be- 
came available the writer this stage the comparison. Their digital 
computer results for the same elasto-plastic problem were then subjected 
the same ductility and factor calculations those used for the authors’ data. 
The results, also plotted Fig. show somewhat better relationship the 


Structural-Dynamic Investigation Fifteen School Buildings Subjected 
Simulated Earthquake Motion,” Blume, Sharpe, and Elsesser, for the State Cali- 
fornia, Dept. Pub. Wks., Div. Architecture (being printed). 


0.2 0.10 0.12 2.20 0.35 0.51 0.48 0.38 
0.3 0.21 0.13 1.62 0.41 0.44 0.52 0.35 
0.4 0.36 0.15 1.42 0.41 0.44 0.53 0.34 
0.5 0.58 0.21 1.36 0.35 0.51 0.52 0.35 
0.6 0.84 0.19 1.23 0.27 0.67 0.37 0.49 
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REF. (6) 


0.20 
REF. (5) 


REF. (6) 


(FIG. 16) 


0.25 0.50 0.75 
NATURAL PERIOD (SEC) 


FIG, 2.—ELASTIC SPECTRA FOR CENTRO 1940 (N-S) 
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curve, thus, providing some substantiation the curve and the Berg- 
Thomaides data. Still another set and perhaps the only other avail- 
able set inelastic data, was subjected the same procedure and plotted 
Fig. The conformity the curve very good. 

Based upon Fig. the following can noted: 


The parameter, period, has little relationship curve agree- 
ment; 

The parameter, damping, seems have little relationship 
curve agreement; 

The compliance this case least, compared the many 
other uncertainties practical earthquake design, suggest that the simple 
reserve energy procedure might applied directly elasto-plastic energy 
design without the need for computer analog analysis. Although the more 
general application complexinelastic yet compared tests 
digital analog results, the writer has, thus far, found good general 
agreement with earthquake history; 

Any elastic acceleration response spectrum might reduced equiva- 
lent inelastic design simply multiplying Rj, determined for any limit- 
ing deflection ratio for elasto-plastic system. (This inverse state- 
ment the third statement.) For example, from Fig. assume the single- 
mass design problem requires Centro 1940-type resistance safety 
factor, the limiting deflection times the yield point deflection, the struc- 
ture’s period 0.5 sec, and the damping critical. The factor 
0.29 regardless damping and the elastic value for Centro 1940, 
damping and period 0.5 sec 0.60g.4 Then, the equivalent code 
design elastic value 1-1/3 times normal stress simply 0.29 0.60g 
0.17g, with energy safety factor. 


noted word caution that there are many other factors 
considered actual design and alsothat these energy procedures for struc- 
tures are recent considerations without the test time. However, the inelastic 
approach research and design seems constitute major assault the 
earthquake problem. Time seems very important also, may attested 
earthquake results Agadir, Arequipa, Mexico City and other places. 

The apparent discrepancy results between the analog and the digital work, 
and, also, between the two elastic system results from the New Zealand ana- 
log, need explanation, course, benefit future research. Some error ex- 
pected all such operations apparent from the complexity major 
earthquake and the sensitivity response. The writer has compared all avail- 
able elastic response data for the Centro 1940 north-south record. These 
efforts several investigators cover range period from 0.10 sec 3.5 
sec, damping from about 20% and, although principally for single-mass 
systems, 10, and mass systems were also included. There was fair 
agreement some areas and considerable variation others. the long 
period, multi-mass range there are onlytwo current sources one which, the 
authors’ electric analog, produces shears less than recent digital source. 
The important point here that the recent higher values, even found 
more accurate, would not change but would instead reinforce the conclusions 
the paper which the multi-mass analog results were originally shown.2 


“Intensity Ground Motion During Strong Earthquakes,” Housner, Of- 
fice Naval Research, Contract N60NR-244, August, 1952. 
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other words, all results data indicate the vast gap between elastic values 
response and code design values even very high assumed damping values. 
The need for realistic energy dissipation design abundantly clear. 

believed that the writer’s staff provided the shadowgraph input data 
for the authors’ analog work although this not certain, since men- 
tioned the authors. Several earthquake records, including Centro 1940 
(N-S), were developed into the form needed for the New Zealand analog and 
furnished the Dominion Laboratory exchange basis, wherein certain 
results were returned. Whether not the Dominion Laboratory has modified 
these graphs not known the writer. 

Comparison previous (1959) from the New Zealand analog with 
the Office Naval Research program results® are shown Fig. and indi- 
cates fairly close agreement both low and high damping values. However, 
subsequent (Fig. 16) analog results, also shown Fig. indicate considerable 
disparity from the previous analog the same shadowgraph were used 
for both series analog tests, would seem that any human errors would 
minor compared other problems. the shadowgraphs were modified the 
Dominion Laboratory, between the two test series, the graph well other 
factors should considered. Two other sets earthquake data were com- 
pared. The Olympia 1949 spectrum, made 1959, the compares 
very well with the record for 20% damping 0.5 sec period, 
which the limit this run. However, the 20%damping spectrum for the 
Centro 1934 (N-S), made the analog, produces accelerations considerably 
greater than that for Housner’s 20% curves® that also 0.5 sec period. 
The Office Naval Research program values, this case, are approximate- 
50% 75% the analog values. 

further question that might considered whether not the digital 
has treated damping beyond the yield point the same manner 
did the authors. The authors apparently combined damping viscosity and energy 
dissipation beyond the yield, which procedure seems logical. 

the writer’s opinion, actual damping values cannot safely assumed 
range from 10% 20% the authors state. Values have often been 
below 10% below for low amplitudes. The use and 20% inthe 
examples this discussion above for comparison purposes only and cover 
range damping. The reserve energy procedure contemplates the 
normal use only reasonable maximum for non-yield and non-damage 
conditions. 

note warning also the sole use the elasto-plastic 
research and design. Many materials and combinations build- 
ing materials are not elasto-plastic even good approximations this 
convenient characteristic. The over-simplification earthquake motion 
simple harmonic years ago was some respects similar over- 
simplifying inelastic response the elasto-plastic. Both assumptions are 


convenient and very useful, but they may vary from actuality. Their limita- 
tions must recognized. 


“Effect Inelastic Action the Response Simple Structures Earthquake Mo- 
Civ. Engrg., 1959. 

“Period Determinations and Other Earthquake Studies Fifteen-story Build- 
ing,” John Blume, Chapter 11, Proceedings World Earthquake Engrg., 
June, 1956, Room 1039, 465 California St., San Francisco, Calif. 
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The discussion herein not intended critical but provide numeri- 
cal data andother data for the apparent differences certain 
analog and computer results. The finding and resolutionof such differences 
itself contribution research and seems proper basis for proceedings 
discussion. The electric analog valuable research tool which has advan- 
tages, (as well disadvantages) over digital and other computers. 

The authors are commended for their pioneering efforts elasto- 
plastic spectral response with the electric analog. hoped that this work 
will greatly extended after further evaluation the procedures and equip- 
ment, and that, eventually, complex inelastic systems can tested well 
elasto-plastic systems. 


ad 
| 
on 
Ts 
| 
‘ter 
A 


DISCUSSION 185 


LATERALLY LOADED THIN FLAT 


Discussion Yervant Terzian 


YERVANT TERZIAN,! ASCE.—The author should complimented 
for concise research paper with important practical application and par- 
ticular for the the results obtained utilizing relatively inexpensive 
equipment. 

Referring the section rectangular plates, the writer carried out se- 
ries tests? 1954, using rudimentary set-up, Fig. witha 
distributed load over the entire surface. The two long sides were always kept 
free, and the short sides were fixed simply supported knife edges. Be- 
cause the difficulties setting the three parallel planes accurately, the 
writer only obtained the curvatures, moments along the axis parallel the 
free long sides. The model used was Plexiglass 1-A, in. in. 1/4 
in., coated with very thin (molecular) coat aluminum vapor vacuum 
chamber. The back the plate was cemented rubber pad with attached 
nylon strings fastened through interconnected levers lead weights, simulat- 
ing uniform load 0.5 psi, all cases. 

The grid consisted series black-nylon wires that were strung from 
lead weights. adjusting the two side plates holding the brass rod, the whole 
set-up was adjusted backwards and forwards and tilted. 

The photographic equipment consisted camera with in. lens 
noted for its perfect spherical abberation correction. 

The lighting arrangement consisted two No. floods with white card- 
boardreflector set in. behind the grid with holein the camera 
lens. Arc lighting did not prove satisfactory. 

Plans and elevations the model mount and the grid are shown Fig. T2. 

For these experiments, small flashlight was used collimate the nega- 
tive and model planes. optical flat was then set the gridplane string- 
ing between two nylon wires. Thus, the grid plane was collimated parallel 
the negative plane. 

Using the presan system, two separate pictures, before load and after load, 

were taken and the grids were measured with filar micrometer microscope. 
the end the experiment, the model plate was cut and milled strips 
in. 1/2 in. order determine the modulus elasticity. The strips 
were cut lengthwise, crosswise, with and without the rubber backing. Load- 
deflection test consisted concentrated load mid-span. The deflection 


October, 1959, William Bradley. 

Assoc., Lawrence Beck and Assocs., Philadelphia, Pa. 

Photographic Method Structural Analysis,” Yerrant Terzian, Master’s 
Thesis, Univ. Pennsylvania, 1954. 
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the center was measured, the results plotted and the modulus elasticity 
determined. The results showed that the rubber backing did not have any ef- 
fect the modulus elasticity. The results obtained were corroborated us- 
ing finite-difference solution, but due the unavailability digital com- 
puter the time, the number unknowns was limited. 

for specific case, the results the Experiment uniformly loaded 
plate, simply supported two edges (short sides) and free two opposite 
edges (long sides), are tabulated Figs. Photographs the negative 
before load and after load enlarged 4.33 times are shown Figs. and 10, 
respectively. The results this series experiments showed that with the 
equipment used, results can ascertained within percentage error margin 


Floodlamp 


Supports 


Surface 


MODEL GRID CAMERA 


FIG, T1.—TEST SET-UP 


10% 15%. commercial semi-automatic devise, involving expensive and 
complicated equipment high accuracy, achieves far better results. 

Mr. Bradley’s work very encouraging from the standpoint deriving 
practical results for complicated problems without resorting this 
cial equipment. could become useful tool consulting office, particu- 
larly with the present tendency towards irregular column positions under 
flat slab. 

would very enlightening the author would delve into the procedures 
followed set the equipment and set the planes the camera negative, 
grid, and model their exact parallel positions. possible source er- 
ror, the suggestion inaccurate collimation the negative, model and grid 
planes may worth discussing. 
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Front View Bection A-A 


Model Mount 


View 
Grid and Tripod 


FIG, T2.—PLANS AND ELEVATIONS, MODEL MOUNT AND GRID SCALE 
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Before Load 


6 in. 

48.31 in. 

50.81 in. 
12x6 


6.38 


0.611 


1.63 except 
supports 


UIRE 
which less than t/3 0.0833 
Check Minimum Radius Curvature: 1/R=0.02022 


=0.02022x0.39 


= 0.00789 
127 in. which more than 96.62 
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After Load 


eee 


Model 


MODEL 
Date March 11, 1954 Loading 0.5 

6.4 0.96 0.04 1.09 

36.45 
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Before Load 


Date March 1954 Loading 0.5 


P= 6 ine 

h= 48.31 in. 


which than t/3 0.25/3 0.0833 in. 
Check Minimum Radius Curvature: 1/R 
0.00726 
137 in. which more than 96.62 
FIG. 


2x6 =6.38 in. 


0.611 


supports 


21.59 1.54 1.06 -0.06 
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After Load 


DIAGRA 


Scales ine Horiz. 1.61 in. Model Vert.= 5.43 in.-lbs. 
Model 


Date March 11, 1954 Loading 0.5 
35.80 0.94 0.06 0.98 
29.56 0.72 0.28 3.15 
0.69 0.31 3.91 5.14 
20.19 0.72 0.28 3.48 
FIG. 


Before Load 


Date March 1954 Loading 0.5 


DESIGN REQUIREMENT 


Check Minimum Radius Curvature:- 1/R=0.02022 
0.00849 
117.75 in. which more than 96.62 


FIG. 
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3 1.64 A= 48.31 ine 
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in. 
& 
1.54 1/J=h¢*p = =0. 
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After Load 


Date March 11, 1954 Loading 0.5 


eee 


MOMENT DIAGRAM 
Scale: in. Horiz. 1.61 in. Model Vert.= 5.43 
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FIG. T9.—BEFORE LOAD 


FIG. T10.—AFTER LOAD 
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With regard large plate deflections, similitude equations would diffi- 
cult satisfy with plastic material corroborating structural steel the 
strain-hardening range. 

would great value obtain influence surfaces electronic digital 
computers using the finite difference solution for rectangular plates with vary- 
ing degrees skew, edge conditions, and continuous spans. 
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BEHAVIOR STRUCTURES DURING 


Discussions John Blume, and Edward Cohen and Samuel Weissman 


JOHN BLUME, author’s paper contains interesting aver 
aged spectral response curves. also brings together, for convenience, sub- 
ject matter that has been treatedin various previous reports. Some 
the subjects have been more extensively considered certainother papers, 
the reference which and related sources might have been interest 
the reader view the general title. 

The average spectrum curves (Figs. and are basedupon the eight com- 
ponents the four strongest ground motions yet recorded. The desirability 
obtaining smooth curves, rather than the jagged curves single spectrum 
for low damping values, recognized. However, sense complacency 
should not developed from the fact that the records used the author for 
this purpose are “the four strongest ground motions yet recorded.” Record- 
ing history only goes back 1933, brief interval for random occurrences 
this nature. would expected with the continued accurate recording de- 
structive ground motion that has been done the Seismological Division 
the U.S. Coast and Geodetic Survey (USGS), that many greater response spec- 
tra will become available with time. The San Francisco 1906 earthquake, 
one example that was not recorded, was earthquake much greater energy 
release and presumably effect structures (or spectral re- 
sponse) than any those recorded thus far. The use averages is, there- 
fore, not yet recommended for ultimate design procedures unless good judg- 
ment allowances are made for the variations intensity, epicentral distance, 
and the geological and soil conditions that would expected with more ade- 
quate sampling earthquake ground motion. 

The variations response for various epicentral distances and soil con- 
ditions has been considered The writer has done considerable 
thinking about this matter and the the greater -intensity close- 
epicenter shocks yet accurately recorded could well create “spike” 
impact effects records that could accentuate short period response con- 
siderably. Moreover, the shocks could have much greater reso- 

nant response long period structures than has been shown spectra thus 
far. The Mexico City and the Taft earthquakes indicated this effect. The au- 
thor’s Fig. is, thus, very interesting the writer who had constructed 


October, 1959, George Housner. 
Pres., John Blume and Assocs. Research Div., San Francisco, Calif. 
Ground Earthquake Motion,” Gutenberg, Bulletin the Seis. 
Soc. Amer., Vol. 47, July, 1957, pp. 221-250. 
“Structural Dynamics Earthquake-Resistant Design,” John Blume, Pro- 
ceedings, ASCE, Vol. 85, No. July, 1958, and Closure, September, 1959. 
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“absolute maximum” spectrum possible spectra, so-called “spectra- 
graph,” the basis response calculations hypothetical impact pulses 
and repeating (resonant) cycles from possible future shocks and the detailed 
study existing spectra, supplemented estimates and plain guesses. The 
purpose was provide limiting value for standard earthquake rating for 
structures all types.3 expected and hoped that such velocity values 
never occur. Fig. shows this “spectragraph” well other spectra in- 
cluding certain the author’s values. Any these can converted ac- 


All curves shown are for critical damping, value that has been sug- 
gested the writer logical value for normal design purposes. Most 
structures, especially buildings, seem have about non-damaging damp- 
ing capacity and more than this cannot, generally, justified the so-called 
elastic range. Much greater energy dissipation capacity may present, but 
this inelastic problem for which other considerations are indicated de- 
pending upon the type motion—steady other and the materials in- 
volved. The damping spectrum has little use design applications. 

The writers “spectragraph” suggested absolute curve for rating purpos- 
the upper line, (a). previously this energy roughly 2.3 
times Centro 1940. Thus, curve (b), 
plotted approximation Centro 1940 (NS). Curve (c), interpolated 
spectral values for and curve (b) are close agreement the 
very short period range which, perhaps, none the curves shown are ac- 
curate. The author’s “average” velocity spectrum the lowest curve and 2.7 
times this his suggested curve for Centro 1940. Although curves (b), (c) 
and the author’s Centro average are very similar shape, there defi- 
nite spread velocity values, much which may due the process 
averaging from records and then factoring back the original; and some is, 
course, due averaging both components for Centro 1940. 

would seem that the writer’s suggested relationship, velocity equals 1.83 
both simple andreasonably accurate for design purposes for Cen- 
tro 1940 (N-S). course, other factors, velocity equals could 
applied for other earthquakes damping values. For values less than 
0.15 0.20 sec, more detailed data will shortly available.5 factor 
applied the author’s average curve Fig. almost perfect coinci- 
dence with the writer’s curve (b) would obtained except the short period 
range where the writer has allowed for the greater close-epicenter shocks than 
have been recorded this manner. 

The curves are Taft 1952, Hollywood 1952, and San Francisco 1957, 
respectively, all which are for damping. From Fig. can seen that 
the short period hump characteristics Taft curve disappear with damping. 
fact, damping, the Taft curve would plot 1.6 times the “average” 
shown Fig. 

The author mentions the effect damping and the value cracking, rub- 
bing, and plastic yield composite building structures and the anomalies be- 


Effects Composite Structures Subjected Earthquake Vibrations,” 
Jacobsen, Dept. Mech. Engrg., Stanford Univ., March 1959. (California 
State Div. Architecture, sponsor). 

Structural Dynamic Investigation Fifteen School Buildings Subjected Simu- 
lated Earthquake Motion,” Blume, Sharpe, and Elsesser, for the State California, 
Dept. Pub. Wks., Div. Architecture. (being printed) 
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tween observed behavior and design values old buildings and discusses 
the damage structure absorbing energy and the survival 
these structures being “not much matter static lateral strength 
question whether the structure can continue absorb appreciable energy 
during its vibrations.” also compares complex, traditional-type building 
simple roof and column structure. These matters have been extensively 
and numerically considered The Alexander Building has been ex- 
tensively studied for many years® and qualitative and numerical values the 
energy dissipation capacity have been 

The author states that the such analysis (overstressed non-linear 
vibration) are not yet form used for design and suggests ap- 
proximate method for assumed elasto-plastic behavior based upon average 
spectral velocity. There general method for aseismic design any type 
inelastic response (including actual failures elements) that form 
for design and also terms that are familiar the structural designer. 
This method involves only simple mathematics and knowledge the force- 
deflection characteristics the structure basic unit same. theforce- 
deflection spring factor characteristics are not known cannot estimat- 
failure, method inelastic design energy dissipation design can 
consistently reliable, matter the analysis. This procedure, 
termed the “reserve energy technique,” has been found simple but ef- 
fective basis for designor predictionof structural response, drift and damage. 
Although far only suggested for use modify analyze structures 
first static codes, seems have interesting possibilities.5,7 
can also used evaluate “plateau” multi-stage resistance, provide 
value against collapse severe emergencies, account for period variations 
with damage, and evaluate condition after damage. 

The authors 18, 19, 21, and are based uponaverage spectral veloci- 
ties. Although the simplicity this often justified, the vari- 
ations that occur velocity with period (Figs. 5,6, the writer’s Fig. 
can appreciable, especially for value that squared obtain kinetic 
energy. would that some nominal damping value, say criti- 
cal, should used for velocity (or acceleration) energy design procedures. 
normal hysteresis loop about both sides the no-load point force- 
deflection curve would account for this, its equivalent, trading energy 
lost the form heat rather than damage permanent set. Assuming 
for this purpose, and using the author’s Fig. values well curve (b), 
damping value. noted that the required energy dissipation values 
vary considerably with period. simple device for improving the relative re- 
sults use the writer’s velocity approximation (Fig. A), modified for the 


“El Centro type” earthquake. Thus, the total energy requirement would 
VT, 0.052 column this basis, the energy dissipation 


the comparisons Table The 0.10 requirement would 


“Period Determinations and Other Earthquake Studies Fifteen-Story Building,” 
John Blume, Chapter 11, Proceedings the World Conf. Earthquake Engrg., 
June, 1956, EERI, Room 1039, 465 California St., San Francisco, Calif. 

Discussion “Electrical Analog for Earthquake Yield Spectra,” John Blume, 

Proceedings, ASCE, Vol. 86, No. June, 1960. 


Fi 
Kg 
if 
4 
eS 
| 
\ 
Bas. 


June, 1960 


200 


i 
x w 

a 

“| 


DISCUSSION 201 


Column represents the energy requirements using the author’s averaged 
Centro 1940 velocities and column the writer’s Centro 1940 N-S ve- 
locities. There variation between the two, but what even more signifi- 
cant, considerable variation either case from constant value such 
0.10 ft-lb energy per pound weight structure. Most buildings, few 
stories, would the 0.5 sec less period range. 

The spectral response curves can valuable design well com- 
paring earthquakes. Code design requirements, although constantly improving, 


TABLE 1.—VELOCITY AND ENERGY REQUIREMENTS FOR CENTRO 1940 


From 


From Curve 


(2) (3) 


(4) 
0.016W 


0.037W 
0.052W 
0.073W 
0.10W 


0.019W 
0.033W 
0.045W 


are not, the writer’s opinion the author’s conclusion, adequate for the 
final design unusual, very slender, very tall, very critical structures. 
Modern aseismic design analysis should deal numerically well intui- 
tively with energy requirements, yield deflection, and damage under various 
degrees severity and probability. 

The research work being done school buildings for the California State 
Division Architecture, indicates that there can considerable variation 
the type and amount earthquake resistance for different buildings designed 
under the same code and, also, that the short period range the spectrum 
(0.05 0.20sec) needs more attention tothe designand analysis the numer- 
ous low buildings. 


EDWARD COHEN,® ASCE and SAMUEL ASCE.—The 
author has presented paper describing the application response 
spectra aseismic design. The proposal have two design criteria seis- 
mically active areas, one avoid damage from frequent earthquakes and the 
second avoid collapse from earthquakes maximum anticipated severity, 

deserves serious consideration. 

Several years ago, the senior writer similar procedure which 
postulated “standard” and “maximum” earthquake. Designs for the “stan- 
dard” earthquake having probable recurrence periodof 

age negligible minor, easily repaired, amounts. For the “maximum” 


Assoc., Ammann and Whitney, Cons. Engrs., New York, 
Structural Designer, Ammann and Whitney, Cons. Engrs., New York, 

“Rational Analysis Earthquake Distortions,” Cohen, Presented the Cen- 
tennial Engrg. Convocation, Symposium the Amer. Geophysical Union, Chicago, 
September 12, 1952. 
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earthquake the structure would allowed incipient collapse. Dis- 
tortions and ultimate strength would the controlling design criteria rather 
than allowable stress. Recent studies and investigations indicate that our know- 
ledge structural strength and behavior under dynamic loads adequate for 
this purpose. 

The “standard” and “maximum” earthquakes would defined synthesized 
ground motions consistent with the geology the vicinity the building 
response spectra for the standard earthquake. However, methods predicting 
the standard and maximum earthquakes appeared beyond the state the 
art. One the implications building codes, which seems condoned 
some extent the author’s discussion, that the maximum anticipated earth- 
quake the maximum experienced earthquake. This contrary general 
engineering practice relative the prediction other recurring physical phe- 
nomena which cannot predicted sound theoretical basis and probably 
unlikely except very active seismic regions. 

Richter has study which shows promise progress 
the direction actual zoning.” Unfortunately, this excellent paper has 
not had wide distribution civil engineering circles, yet. Richter reports 
that proposals similar that the senior writer were made Russia 
high intensity regions (greater than Mercalli VIII) being divided into areas with 
probable recurrence times centuries, decades, years. Although 
such “risk zoning” still some time off, Richter has prepared regional map 
(Fig. probable maximum intensity for the United States which should 
considerable interest the profession. 

may ofinterest note that the ground motions associated with nuclear 
weapons bursts are currently being considered the design protected fa- 
cilities. The problems encountered are, many ways, similar those 
earthquake engineering and procedures the type suggested the author are 
currently being used analyze structures and shock isolation systems. How- 
ever, the physical characteristics the ground shocks differ from those 
earthquake. 

Ground motionfrom nuclear bursts may induced twodistinct process- 
es: (1) air-induced shock, and (2) ground-transmitted shock. The ground- 
transmitted shock direct hydrodynamic energy into mechani- 
cal energy the vicinity the explosion. Nuclear-weapons-test data has in- 
dicated that the air-induced shock effects are substantially larger than the 
ground-transmitted effects for most current design pressure regions. 
The air-induced shock caused the blast wave traveling over the ground 
surface. The characteristics the blast wave terms intensity, propa- 
gation velocity, and variation withtime are functions the weapon size, height 
burst and distance from ground zero. This air blast wave impulse 
the ground surface results predominantly single pulse ground 

Fig. shows records typical vertical versus time 
gravity units for the and depth below the ground surface. These 


“Seismic Regionalization,” Richter, Calif. Inst. Tech., Contribution No. 
897, December, 1958. 

“Underground Shock Environment Data and Application the Design Under- 
ground Structures,” Cohen and Weissman, Proceedings the 28th Symposium 
Shock, Vibration, and Associated Environments, Office the Dir. Defense, Re- 
search and Engrg., Washington 25, C., February, 1960. 
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were obtained during nuclear the Nevada Test for 
kiloton weapon the 200 psi over-pressure range. seen Fig. the ac- 
celeration records are characterized single sharp peaks acceleration 
(182 and 54.2 gravities) preceded and followed minor disturbances. The 
duration the peak downward acceleration pulse the order millisec 
for the 1-ft depth and millisec for the 5-ft depth. The total durations are 
about 0.18 sec which are approximately equal the air blast duration. These 
curves are free-field accelerations both the vertical and horizontal 
directions that there rapid attenuation acceleration and decrease 
frequency with depth. The single sharp acceleration pulse corresponds 
ground displacement which also essentially predominant single displace- 
ment pulse. This differs from typical earthquake motion that the frequency 


Air Blast Arrival 


Depth 


0.200 0.300 
Depth 


TIME (Seconds) 


FIG, 2.—VERTICAL FREE—FIELD GROUND ACCELERATION VERSUS 
TIME, NEVADA TEST SITE; W=40 KT, p=219 PSI 


and intensity acceleration much higher, having single sharp peak, rapid 
decay and considerable attenuation with depth the upper layers. 

Fig. shows vertical displacement shock spectra recorded the Nevada 
Test for kiloton weapon yields 116 146 psi pressure re- 
gions, multi-frequency reed gages with approximately one-half per cent 
critical damping. Peak displacements the gages were recorded near the 
ground surface (2-ft depth) the free field and shallow buried structure. 


“Ground Acceleration, Stress and Strain High Incident Overpressures,” 
Swift, and Sauer, Project 1.4, Operation Plumbbob, ITR 1404, Stan- 
ford Research Inst., Menlo Park, Calif., 1957, (Classified). 

“Ground Shock Environment—Measurement and Application,” Barton, GM- 


TR-293-A, Space Tech. Lab., The Ramo-Wooldridge Corp., Los Angeles, Calif., De- 
cember 11, 1957. 
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indicated the spectra there considerable attenuation the free-field 
responses the shelter. Fig. vertical acceleration spectrum derived 
from average displacement responses Fig. The responses increase 
the frequency increase reaching peak about 150 cps and decreasing the 
higher frequency range. significant note that the range most earth- 
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FIG, 3.—VERTICAL DISPLACEMENT SHOCK SPECTRA, 


quake spectra would fit into the small area indicated near the originof the ac- 
celeration spectrum curve. Although data limited appears that response 
records particular ground shock gage not exhibit complicated form 
the earthquake spectra. This the result oscillators all frequencies 
responding the predominant single pulse, whereas earthquake spectra may 
representative resonant responses the particular frequency range close 
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Frequency, cps. 
FIG, 4.—AVERAGE VERTICAL FREE—FIELD GROUND ACCELERATION 


SPECTRUM, NEVADA TEST SITE; W=10TO KT,p 116 PSI 
146 PSI 
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Soil Undisturbed Region 
(Free-Field) 
FIG, 5.—SCHEMATIC REPRESENTATION STRUCTURE 
GROUND SHOCK EARTHQUAKE FIELD 
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the frequencies the ground vibration. For the same reason, structure 
damping not important this case the case earthquakes. 

important note that spectra obtained from free-field ground motion 
data (motion not affected the presence structure) may not the actual 
motion the ground the base multi-story building surrounding 
buried structure.15 This particularly true for the high frequency portion 
the spectra that the interaction the structure and the ground will, ef- 
fect, modify the actual vibrations induced the ground. blast induced shock 
this effect magnified the sharp variation the motion with depth. Thus, 
ground motion data obtainedin the basement particular structure may not 


applicable structures different mass, configurations and frequency 
even the same area. 
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FIG, 6.—DISPLACEMENT AND VELOCITY PULSE 


Fig. 5is schematic representation structurein general ground shock 
earthquake field. The motion the soil mass adjacent the structure will 
modified from that the free-field the interaction the structure and 
the soil. Although the free-field response spectrum which exists the undis- 
turbed region cannot used directly for analysis the structure, may 
applied the overall system, including the adjacent soil mass the systemis 
elastic. used for elastic structures where the ground 
interaction may neglected approximated. For elasto-plastic, non-linear 
discontinuous structures free-field spectrum, without further analysis, 
has only limited value. synthesized ground motion curve consistent 
with both the spectrum and observed motions useful. 


“Comprehensive Numerical Method for the Analysis Earthquake Resistant 


Structures,” Whitney, Anderson, and Salvadori, Proceedings, ACI, 
September, 1951, Vol. 48, pp. 5-28. 
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Such ground motions have been synthesized the case nuclear ground 
shock approximating the ground motionas single velocity 
Fig. shows displacement and velocity pulse which was derived 
from the average response spectrum Fig. This pulse, applied series 
single degree freedom oscillators, would result the peak responses 
the spectrum for the frequencies below 100 cps which includes the range 
application most structures. velocity pulse having the peak 
value may derived conform more closely the higher frequency area 
required. The displacement curve the same type observed 
Using such pulse consistent analysis may performed numerical inte- 
gration taking into account elasto-plastic, non-linear and discontinuous be- 

The longer duration and more periodic nature the earthquake ground mo- 
tion prevents the use single pulse and complete time history the ground 
motion would have considered inan analysis this type. ground 
displacements varying frequency could determined from typical record- 
earthquake motion and spectra. 

discussed and illustrated the author suitably designed structure 
which stressed beyond the yield point dynamic loads will absorb ap- 
preciable amount energy without collapse. also appears that because 
the increased earthquake intensity required ductile structure than 
damage it, the strength requirements the earthquake may not 
greatly exceed may even less than those the “standard” earthquake. 
Studies blast have indicated that feasible substantially 
increase the blast overpressure required produce collapse the conven- 
tional structures minor modifications the initial design fully mobilize 
the inherent lateral strength and ductility the construction. However, 
the validity design based upon dynamic analysis function not only 
the accuracy the calculations but also the consistency the actual mass, 
stiffness, and strength magnitudes and distributions with those assumed. Where 
large plastic distortions are allowed unanticipated increases variations 
these quantities may catastrophic deficiencies. 


Developed Chobotov, Space Tech. Labs., Inc., Also, See Reference 

“Impulsive Motion Elasto-Plastic Shear Buildings,” Cohen, Levy, and 
Smollen, Transactions, ASCE, Vol. 122, 1957, 293. 

“Design Balst Resistant Construction for Atomic Explosions,” C.S. Whitney, 
Anderson, and Cohen, Proceedings, ACI, Vol. 51, March, 1955, pp. 589-683. 
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